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ABSTRACT 


A subsonic wind tunnel investigation was conducted at the Naval Postgraduate 
School Wind Tunnel Test Facility to examine the effects of grid generated turbulence 
on the high angle of attack aerodvnamics of a vertically launched missile. Four turbu- 
lence generating grids produced turbulence with length scale to missile model diameter 
ratios of 1.05, 0.89. 0.62, and 0.15 with respective turbulence intensities of 3.3195, 
2759. 1.9860, апа 0.47%. The test model is representative of current cruciform wing 
nussiles with low aspect ratio wings and long root chords. The tangent ogive nose has 
a nose fineness ratio of 2.29. Three separate body configurations were tested with and 
Without the turbulence generating grids. One configuration was a bodv in isolation and 
the other two were wing and tail configurations at 0° and 45° roll angles. All test runs 
were conducted at Re,=1.1x1@ over an angle of attack range of -5° to 95°. Results in- 
dicate that as turbulence length scales approach body diameter size. the angle of attack 
for the onset of asymmetric vortices was delayed and the side force magnitude was re- 
duced. The vortices generated bv the nose of the missile continue to dominate the 
afterbodv vortices for body configurations with and without wings regardless of the 


turbulence conditions. 


THESIS DISCLAIMER 


Ihe reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of computational and logic errors, 
they cannot be considered validated. Any application of these programs without addi- 
tional verification is at the risk of the user. 
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A. BACKGROUND 

Тһе development of missiles capable of operating at high angles of attack has 
prompted numerous investigations and experiments in high angle of attack research. 
Much interest is due to the demands for highly maneuverable missiles at supersonic 
speeds. In cases where the missile is launched from the ground or at sea from a ship. the 
missile’s initial velocity for all practical purposes is zero. From this initial state, the 
missile must flv over a wide range of Mach numbers from subsonic to supersonic speeds. 
The aerodvnamic characteristics of the missile may vary greatlv as it transitions through 
these operating speeds. All surface-to-air missiles follow a three phase operation: 
launch. midcourse. and terminal. Each phase represents a different missile flight regime 
owing to the changing flight control requirements and airflow about the missile. thus 
presenting a challenge to the mussile designer who must develop a missile capable of 
Operating throughout each regiine. 

Vertical launch for missiles is by no means a novel idea. It has been the primary 
means of launching large rockets and missiles from land. Vertical launch of surface-to- 
air missiles at sea is a recent development. Conventional launch systems utilize a 
trainable launcher consisting of a rail or tube adjustable 1η azimuth and elevation. The 
vertical launch surface-to-air missile (VLSAM) uses a canister container which stores 
and launches the nussile. The canister system eliminates the cumbersome trainable 
launcher and associated specialized storage facilities and loading mechanisms. Other 
advantages of the VLSAM are the 360? coverage provided without interference from the 
ship's superstructure. a higher firing rate since launcher availability is no longer a factor, 
and space and weight savings [Ref. 1]. 

VLSAMs have their own unique trajectories, allowing them to point at their targets 
after launch. Two flight control methods are currently used: thrust vectored control, and 
aerodynamic surface control [Ref. 2]. Thrust vectored control 1s used with short range 
missiles which must pitch over quickly to engage targets at close range. Long range mis- 
siles maneuver with aerodynamic control surfaces. In each case, as the VLSAM exits the 
canister it enters the open ocean environment at low velocity and 1s subject to potentially 
significant crosswinds. The resultant of the missile and crosswind velocities 1s a potential 


high angle of attack flow about the missile [Ref. 3]. This high angle of attack flow may 


cause asymmetric vortices to form on the missile nose and afterbody which induces side 
forces as discussed later in this chapter. These side forces could give flight control 
problems to the VLSAM at launch. | 

The launch environment may contain some degree of turbulence, both from the at- 
mospheric boundary laver, and the airflow about the ship. Understanding the effect of 
this turbulence on the aerodynamic characteristics of a VLSAM during launch is the 
goal of a continuing effort of research and experiments conducted at the Naval Post- 
graduate School (NPS). In earlier work, Roane [Ref. 3] developed a system of modelling 
the flowfield turbulence by a series of four grids to generate turbulence in the NPS low 
speed wind tunnel. The goal of this thesis is to experimentally examine the effects of this 


turbulence on the asymmetric vortex induced forces generated on a VLSAM model. 


B. HIGH ANGLE OF ATTACK AERODYNAMICS 

Flow separation from the body, wing, and tail surfaces is an important phenomenon 
encountered in high angle of attack aerodynamics since asymmetric vortices may form 
to induce side forces. Primary factors which influence flow separation from the nussile 
body are: the angle of attack, crossflow Reynolds number (Re,) !, nose geometry. and 
nose length to diameter ratio (also called the nose fineness ratio). Secondary factors 
include roll angle and roll rate, free stream turbulence, surface roughness, acoustic en- 
vironment (pressure waves), and vibration [Refs. 4, 5, 6]. This section discusses how 
the factors mentioned above influence the mechanisms responsible for vortex generation 
and disposition. 

l. Vortex Regions. 

As à slender unvawed bodv changes angle of attack « from 0° to 90°. the airflow 
over the missile body transitions through four regimes as described by Ericsson and 
Reding [Ref. 7: pp. 246-247]. These regimes [Figure 1] are characterized by the chang: 
ing vortex patterns and the resulting induced in-plane (normal). and out-of-plane (side) 
forces and moments which act on the missile at zero side-slip and at high angles of at- 
таск 21Кеі5, 2, 4, 6, 9 |. λαο ο ου. 

бу - the angle of attack at which steady svmmetric vortices are formed. 
«ду - the angle of attack where steady asymmetric vortices are formed. 


осу - the angle of attack where unsteady vortices are formed. 


1 The crossflow Reynolds number is set by the base diameter of the missile. 


2 Without wings and tails. 
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Figure 1. Vortex Generation Regimes: (Ref. 10] 


a. Regime I - [0? <a < as] 

The axial flow dominates and flow around the body remains attached for 

pointed noses. If the nose is blunted. a separation bubble may occur (discussed below). 
b. Regime Il - (>: Sa <a] 

The crossflow begins to push the boundary layer to the leeward side where 
it separates and a system of two or more pairs of symmetric vortices are formed. These 
vortices originate at the nose of the body. For long bodies, more pairs of symmetric 
vortices are generated along the body length. The number and strength of the symmetric 
vortices increase with angle of attack. Symmetric vortices typically start at o..35° 
| 4) 

c. Regime III - {Zy S a < zu] 

Crossflow begins to dominate, shedding asymmetric vortices which induce 
out-of-plane or side forces on the body. These asymmetric vortices are relatively steady, 
but may change from side to side at higher angles of attack near the maximum side force 
magnitude [Ref. 11]. In some cases the side forces may exceed the normal force. The 
vortices which are the most asymmetric yield the highest side force magnitude 
[Ret [5], i his regime tipicas starts at oa mO Req] 

d. Regime IV - [wc < х < 90°] 

The crossflow dominates completely and flow separation becomes unsteady 
and dependent on Revnolds number. At the critical Reynolds number the boundary 
laver separates then re-attaches as a turbulent vortex sheet (discussed later) Ihe 
boundarv laver 1s shed as a random wake at supercritical Revnolds numbers. or a von 
Karman vortex sheet at subcritical Reynolds numbers. Normal force magnitude will 
plateau and side force magnitude decreases to zero. This regime has been observed to 
begin at a.a 60? [Ref. 4]. 

Asymmetric vortex induced side forces are dependent on the strength and 
disposition of the vortex patterns which are, according to Ericsson and Reding [Ref. 7], 
a result of boundary layer transition on the surfaces of the missile. The analogy most 
often used to describe flow separation deals with boundary laver transition on a two- 
dimensional cylinder in a normal crossflow. Three dimensional effects due to nose and 
afterbody geometrv must also be understood as they dictate the formation of the vortex 


structure: 





2. Two Dimensional Crossflow. 


Airflow over the nussile body can be divided into normal and axial components. 
The axial fiow component follows along the missile bodv length and the crossflow is 
essentially a two dimensional flow normal to a cylinder. The crossflow analogy provides 
information for cylinder lift and drag which act in the direction of the crossflow. Side 
forces may exist at right angles to the crossflow depending on the type of flow separation 
on opposite sides of the cylinder. Mechanisms behind boundary laver transition and 
separation provide an explanation for flow separation and subsequent asymmetric 
vortex generation effects. 
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Figure 2. Two dimensional crossflow about a cylinder: [Refs. 13, 14]. 


The crossflow Reynolds number is the primary factor which influences the sep- 


aration point of the boundary layer. Viscous flow, surface roughness, and turbulence 


are other factors which influence boundarv laver separation. Achenbach classined the 
flow around a cvlinder into four distinct regions deliniated by differing tlow separauon 


and drag behavior [Ref. 13]. These regions are shown in Figure 2and discussed beiow. 
In the subcritical range, the boundary layer 1s laminar and flow separation oc- 


curs close to to the lateral meridian where @. defined as the angle from the direction of 
the crossflow, varies from 80° to 90°. Asymmetric subcritical separation occurs at 
@=90°, where the side force just begins to reach a noticeable magnitude. 
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Side force to normal force ratio: [Ref. 7: p. 260]. 


Figure 3. 
When the Reynolds number enters the critical range, the laminar boundary laver 


separates from the body at @%90° and re-attaches as a turbulent boundary laver which 
separates again from the surface at @=140° . Between the laminar separation and the 


turbulent re-attachment, a laminar bubble is formed. In this region, the flow separation 


can easily change back and forth from critical to subcritical for small changes in 
Reynolds number. The changes occur on opposite sides of the cylindrical section and 
when critical separation exists on one side of the body and subcritical separation on the 
opposite side. a large difference in ó from opposite transition points is possible. When 
this critical and subcritical combination occurs. an asymmetry is observed in the gener- 
ated vortices. The pressure difference on opposite sides of the lateral meridian will 
produce the side force. For the largest difference in @. the vortices will be at maximum 
asymmetry, and the side force at the highest magnitude as seen in Figure 3. The 
asvmmetric vortices mav also become susceptible to altering back and forth, hence the 
change in the direcuon of the induced side forces. Maximum side force occurs at the 
critical Reynolds number. making it an important parameter. [Ref. 15] 

Next. as the Revnolds number increases into the fully critical range. laminar 
transition moves towards the direction of the crossflow. Turbulent separation will move 
forward to 6 < 140° and the laminar bubble is no longer formed. Asynunetric critical 
separation occurs simultaneously on opposite sides of the body. The asymmetric vortices 
produce nearlv equal surface pressures at the lateral meridian because the separation 
asymmetry only affects the surface pressures at ¢@ > 140°. In this range the vortices are 
ineffective at producing a significant side force and a sudden decrease in magnitude oc- 
curs as shown bv the curve in Figure 3. The boundary laver thickens and viscous 
properties of the fluid are the dominant factor. Drag reaches a minimum. 

Finallv, as the Revnolds number increases further to supercritical and into the 
transcritical range. the laminar transition point works towards @=0° and turbulent sep- 
aration occurs at @= 100°. The asymmetric transcniucal separation point moves towards 
the lateral meridian where the asymmetric vortices once again produce a significant side 
force. Drag vill continue to rise and settle out at a level lower than subcritical condi- 
uons. [Refs. 7: pp. 247-261. 13] 

The Revnolds number dictates the greatest influence on the normal force and 
drag characteristics, especially within the critical range. Maximum normal forces on the 
missile should occur when the boundary layer flow is completely laminar or completely 
turbulent, because crossflow drag is high. Surface roughness and viscous effects are also 
present. Experiments by Lamont [Ref. 16] and Wardlaw [Ref. 17] show that maximum 
asymmetric loads follow the same behavior as the normal force, with the lowest vortex 
induced side forces occurring at transcritical conditions. The largest lift to drag ratio 


occurs at the critical Revnolds number where drag is a minimum. Maximum vortex 


asymmetry has been observed to occur at the critical/subceritical Reynolds number ex- 
plained above. Thus the maximum side force to normal force ratio may be expected to 
be found at the critical Revnolds number. [Ref. 14] 

3. Three Dimensional Vortices. 

The initial separation point of the shed vortices occurs at the nose of the slender 
body. Many researchers have observed that these vortices shed from the nose will dom- 
inate other vortices shed along the body length [Refs. 18, 19, 20]. Thus, the nose ge- 
ometry becomes an important factor in vortex generation and disposition. Axisymmetric 
missile nose geometry. falls into two shapes: cones and ogives, both pointed and blunt. 
Noses are dimensionalized for comparison by their length to base diameter ratio or 
limmeness- ratios (1 g map 

For pointed conical and ogive noses, observations indicate that a4, occurs as a 
function of the semi-vertex angle 0, . At all Mach numbers, asymmetric vortex shedding 
begins when the the angle of incidence o 1s greater than the apex angle of the nose, where 


0. = 0c for a conical nose. 
AVR 20A (1) 


For a tangent ogive nose3, with length 4; and base diameter d:- 


N 
5 / 
OM — (2) 
ΠΠ. 
7 0 25 
A close approximation for slender bodies is: 
c πα (3) 
Ix 


Blunting the nose may cause a second type of asymmetric vortex shedding to 
develop at the rear of the nose. The blunt nose causes a separation bubble to form on 
the leeward side of the nose. Asymmetric vortices originate at the rear of the blunted 
nose, and tend to dominate the induced effects from the afterbody vortices [Ref. 21]. 


The asvmmetric vortices generated by a blunted nose do not alternate from side to side 


3 A tangent ogive nose is constructed from a constant radius arc with the center of the arc lying 
in the plane of the base of the nose and some distance away from the nose axis. 


as readily as those generated a pointed nose. Side force magnitudes are also reduced. 
Afterbodv vortices from a blunt nose missile remain closer to the body and are easilv 
shed throughout the missile length [Ref. 4]. Ihus, blunting the nose will reduce the ef- 
fects of vortex asvmunetry and the induced side force. 

Nose fineness ratio has been shown to have an effect on the asymmetric vortex 
induced side force [Refs. 18, 19, 22]. In general, side force magnitude increases with 
larger nose fineness ratio [Ref. 18]. As the nose fineness ratio increases, the nose apex 
angle decreases and the angle of attack for the onset of asymmetric vortices will de- 
crease. Thus, a missile with high nose fineness ratio may become more susceptible to 
induced side forces at a lower angle of attack than a missile having a lower nose fineness 
ratio. Jorgensen demonstrated that decreasing the nose fineness ratio 1s more beneficial 
In reducing side force than blunting the nose [Ref. 1$]. 

Тһе presence of the afterbodv at a constant roll angle has almost no effect on 
vortex generation from the pointed nose of slender bodies [Ref. 23]. The first pair of 
vortices generated bv a pointed nose will set the pattern for vortices shed along the bodv 
[Ref. 24]. The afterbody will also shed asymmetric vortices, but these will not effect the 
nose vortices generation, although the afterbodv forces will contribute to the overall 
magnitude of induced forces from the nose-bodv combination. 

Should the afterbodv roll angle change. with nose roll angle remaining constant, 
the afterbodv vortices will change the side force magnitude and direction as shown in 
research by Kruse [Ref. 25]. Since the nose and body of a missile do not rotate inde- 
pendently of each other. the vortex system generated by the nose should still dominate 
afterbody vortices on a axisymmetric missile configuration. 

Experimental results for pointed ogive-cvlinders with and without afterbodies 
show that small variations in roll angle alter the asymmetric vortex structure and side 
forces. Researchers believe that surface imperfections and deviations in the nose 
axisymmetric geometry are the sources of this phenomenon. The direction and magni- 
tude of the side force becomes unpredictable, a detriment to anyone attempting to gen- 
erate a model. [Ref. 26] 

The preceding section on the two dimensional flow separation analogy discussed 
four flow separation mechanisms. Keener identified that three of these flow separation 
mechanisms occurred simultaneously along the length of a nose and forebody at angles 
of attack greater than 20° and transcritical Reynolds number of Re, =0.8x10° . These 


mechanisms were observed from surface flow patterns on tangent ogive noses. 


]. Primary laminar separation near the nose. 


το 


Primary transitional separation at the nose midsection consisting of a lamunar 
bubble. 


3. Primary turbulent separation near the base of the nose. 


The extent of the transitional separation region grew with increasing angles of attack 
and constant Reynolds number while the turbulent separation region decreased. 
[Ref. 11] 

Ongoing investigations have not fully revealed the exact cause of vortex 
asymmetry. The general consensus among high angle of attack researchers 15 that vortex 
asymmetry is induced by the location of boundary layer separation [Refs. 7, 9]. Keener 
also considers the possibility that vortex asymmetry may be caused by an inviscid in- 
stability in the pair of separated vortices. [Ref. 11] There are many other three dimen- 
sional effects which must be considered to complete a flowfield analysis which include 


roll, pitch, and yaw rate; and vaw angle. 


C. LIFTING SURFACE EFFECTS 
1. Wings. 

Compared to aircraft, missiles generally use low aspect ratio wings. Mutual ef- 
fects from a wing-body combination are important since in some cases the mussile wing 
span approaches body diameter. 

As discussed earlier, nose vortices dictate flow behavior over the missile body 
at high angles of attack. In turn, the effects may also be felt by the wings. Without 
wings, the nose and body vortices tend to move away from the body with increasing 
angle of attack. The addition of wings causes the vortices to take a different trajectory, 
moving them closer to the body. As the vortices move closer to the body, the result 1s 
comparable to increasing the effective angle of attack causing unsteady asymmetric 
vortuces. || Relea 

A significant fraction of the normal force from a missile wing 1s proportional to 
the square of the angle of attack. This is attributed to the vortex lift effects from flow 
around a sharp leading edge delta wing 4. As vortices curl around the wing leading edge, 
a low pressure area forms on top of the wing providing lift as shown in Figure 4. Since 
a cruciform missile has an equilateral wing arrangement, equal lift may be provided in 


the normal and side directions depending on the wing angle of attack to the flow. For 


4 A common wing planform shared by numerous missiles. 
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wings with low aspect ratio, a large proportion of the lift generated bv the wing vill be 
gdueuo vendes lit. IRef. 28] 
The net eílect of the wing-bodv combination seems to be à reduction in the ef- 


[ecuve angle of attack for the onset of asymmetric vortices and side forces. 





Figure 4. Flow over a delta wing: The left drawing shows the vortices curling 
over the sharp leading edge for low Revnolds numbers. The right draw- 


ing is the cross section B-B. [Ref. 29]. 


2. Strakes 

For an aircraft, the addition of strakes is to improve the aerodynamic charac- 
teristics by enhancing lift at high angles of attack. Vortices generated bv flow separauon 
from the aircraft's forebody interacts with the flow over the main wing. The addition of 
strakes will generate additional vortices improving the interaction between bodv vortices 
and airflow over the wings. The effects of wing strakes arc different on a low aspect ratio 
wing missile since nussile forebodies are comparatively longer and wing half span 1s 
smaller when compared to body radius. The flow becomes separated from the body at 


a moderate angle of attack before passing over the strakes [Ref. 30]. The strakes ona 
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missile may lessen the effect of asvmmetric body vortices since the strakes would cause 
interference with the crossflow component around the body. The real benefit of incor- 
porating strakes with low aspect ratio wings 1s from the strong vortices produced by the 
strake which enhance vortex lift effects. 

3. Tails 

The tails of a missile have little influence on flowfields on the forebody, espe- 
cially at low angles of attack. At high angles of attack interference on the flowfield 
caused by wings and the afterbody mav effect tailflow which depends a great deal on the 
placement of the wings, the length of the afterbody between the wings and tails, and the 
missile's angle of attack. For higher angles of attack, nose and wing vortices may be 
sufficiently displaced away from the body to have onlv a slight effect on the tailflow. 
Should a large distance between wings and tails exist, afterbody vortices may donunate 
the airflow over the tail surfaces. Hence, most nussile configurations with long 
afterbodies utilize canard or wing control while missiles employing tail control have 
short afterbodies and low aspect ratio wings with long root chords. The additional 
asvmmetric vortices from the afterbodv will also contribute to the overall side force on 
һе пие е | 


D. TURBULENCE 
1. Intensity and Length Scale 

Turbulence denotes an irregular fluctuation or disturbance of a fluid changing 
motion when superimposed on a flowfield. In terms of scale (size and energy), the 
flowfield mav be magnitudes higher than turbulence. The disturbance exists when fluid 
particles overcome the constraining effects of viscosity and are able to fluctuate in the 
flowheld. Turbulence exists for large Revnolds numbers where longer flowfield length 
scales are present. A horizontal wind or crosswind may be considered as a large scale 
flowfield and turbulence becomes a smaller component of the crosswind. Further, when 
a flowfield contains turbulence, it can be considered a turbulent flowfield and so а 
crosswind can be though of as a turbulent crosswind. It should be understood that anv 
wind for practical purposes is turbulent, but since crosswind size and energy scales are 
magnitudes larger than the turbulence scales contained in the crosswind, we can treat the 
crosswind as homogeneous and two dimensional. 

As a large scale homogeneous flowfield, the crosswind can be viewed as exerting 
a uniform flow over the missile body. Flow separation mechanisms discussed earlier for 


a two dimensional cvlinder in a steady flowfield apply in the crosswind. Understanding 


the effects of small scale turbulence on the flow separation mechanisms is important 
since vortex generation and asymmetry determine side force direction and magnitude. 
Small scale turbulence effects are dependent on turbulence length scale and intensity 
which will be discussed below. Turbulence intensity and length scale components exists 
in the streamwise and spanwise directions relative to the flowfield. For purposes of this 
and following discussions onlv the streamwise component in the flowfield will be dis- 
cussed. 

Turbulence intensity is the measure of the relative magnitude of velocity fluc- 
tuations in the flowfield. Mathematically, for a horizontal flowfield or crosswind, the 
turbulence intensity T, is the ratio of the root-mean-square (rms) velocity fluctuation 


u , to the mean velocity component in the flowfield, 0... 





T,2— (4) 


Since the turbulence intensity measures the magnitude of velocity fluctuation, a higher 
Intensity indicates higher kinetic energy and correlates to more turbulent flow. 

Turbulence length scale L, describes the time average measure of the fluid dis- 
turbance eddies. The measure is time averaged since these disturbances constantly 
change size. The size of the disturbances are determined by external conditions associ- 
ated with the flow. An example of large scale turbulence are flow disturbances caused 
bv flowfields over large buildings in a city, while small scale turbulence is generated by 
Memecnrouen mesh screens. [Ref. 31: p. 557] 

Turbulence length scale to body size ratio may dictate the manner in which 
turbulence affects the VLSAM flowfield. The turbulence length scale can be compared 
uoumssuesleueth. L.: EL, , or missile diameter, L: Li- When the turbulence length scale 
Is many times larger than the missile length, L,>L,, the effect on the missile is somewhat 
like a uniform time independent steady-state flowfield, especially at low turbulence in- 
tensities. The effects of the flowfield on vortex development are largely dictated by those 
conditions and factors discussed earlier for a two dimensional cylinder. When the tur- 
bulence length scale 1s comparable to the body length, L,zL, , the flowfield ts distinctly 
non-uniform and may cause unwanted rolling, pitching and yawing motion [Ref. 32]. 
In both cases, the missile body motion results from the actual flowheld conditions about 


the missile body. 
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Perhaps the most significant turbulence condition is when the turbulence length 
scale is smaller than the body. L,«L, , and most important, when it is comparable or 
smaller than missile diameter, L,ZL,. The small scale turbulence present in the flowfield 
mav change the boundary lavers present on the nussile surface, particularly when the 
turbulence length scale is the same size or smaller than the boundary laver thickness. 
An important parameter for small scale turbulence is the turbulence intensity, since the 
energv in the turbulence mav trigger the boundary laver. In contrast to larger turbu- 
lence scales where the primary influence is on the flowfield, the effects of small scale 
turbulence 1s felt on the boundary laver and flow separation on the surface of the missile. 
RET] 

Both large and small scale turbulence length scales actually coexist in the 
flowfield since the turbulence length scale is the mean of the disturbances in the 
flowfield. The combination of different length scales mav add complexity to the qual- 
itative analvsis of the flowfield. The cascade effect explains how turbulence in varying 
length scales exists in the flowfield. As the disturbances in the flowfield experience 
strain, the disturbances will break up into finer disturbances of smaller scale and de- 
creasing energy. The process repeats itself or cascades until the small disturbances 
eventually disappear owing to viscositv. When the disturbances decrease 1n size and en- 
ergy, their individual intensities will decrease at a faster rate. Thus, the turbulence length 
scale average will be biased towards the larger turbulence scales in the flowfield and gives 
an accurate quantitative representation of all the length scales within the flowfield since 
the energv and intensitv of the larger scale turbulence predonunate. [Ref. 34] 

2. Boundary Layer Effects 

Turbulence length scales on the order of boundary layer thickness may change 
the boundarv laver on the missile surface. This small scale turbulence will be referred to 
as boundary layer scale turbulence, and although three dimensional and inhomogeneous, 
can be visualized as two dimensional crossflow components over a cylinder. When the 
boundary layer scale turbulence in the crossflow reaches the missile surface boundary 
laver, it contributes energy as it dissipates into the boundary laver.5 This turbulence 
interaction results in an earlier laminar transition and delays separation of the turbulent 
boundary laver. This causes flow separation to occur below the critical Reynolds num- 
ber for laminar flow, postponing asvmmetric vortex shedding. Revnolds number effects 


dominate at low turbulence intensities for smooth surfaces. When surface roughness in- 


5 In accordance with the laws of momentum and energy conservation. 
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creases, the boundary layer separation becomes less dependent on Reynolds number and 
more dependent on turbulence intensity. The net effect postpones the side force induced 
by the asymmetric vortices and allows the nussile to operate at higher angles of attack 
without induced side forces. [Refs. 33, 35] 

The application of boundary layer scale turbulence to add energy to the 
boundarv laver is a well known method of inducing earlier flow transition and delayed 
flow separation to eliminate or reduce the effects of asymmetric vortex induced side 
forces. Among the numerous methods implemented to add turbulence to the boundary 
laver are placing roughness strips in Various patterns at noses of the missile models uti- 
lized in wind tunnel experiments. [Refs. 11. 26] 

3. Vortex effects 

When the vortices from the missile body are formed. they become detached from 
the bodv at high angles of attack. The interaction of turbulence on the order of the 
vortex size may give additional insight to the disposition of asvmmetric vortices. This 
small scale turbulence will be referred to as vortex scale turbulence and is magnitudes 
larger than boundary scale turbulence. Vortex scale turbulence may be compared to 
body diameter size L,zL,. Objects in the flowfield, much larger then roughness strips, 
тау бе used to generate vortex scale turbulence. Vortex scale turbulence mav be gen- 
erated by strakes. support hardpoints, or inlets. The effects of vortex scale turbulence 
may influence the position of the asymmetric vortices and may contribute dissipated 
turbulence energies as with boundary laver turbulence. Since the relative positioning 
and magnitude of the vortex structure governs the pressure distribution on the surface 
of the missile. the vortex scale turbulence could dictate the position and magnitude of 
the asymmetric induced side forces. Research into this phenomenon has been con- 
ducted as another means of alleviating asymmetric vortex induced side forces. Methods 
include adding nose strakes, nose helical vortex generators, and nose booms to nussile 
models. [Refs. 11, 18, 36] 


E. VSLAM LAUNCH ENVIRONMENT 
1. Ihe Marine Environment 
Conditions in the marine atmospheric environment may have a significant effect 
on the VLSAM. In the atmosphere, turbulence conditions exist almost everywhere, in 
particular, within the atmospheric boundary layer (ABL) which is the lowest part of the 
atmosphere [Ref. 37: p. 179]. The conditions in the ABL are the result of atmospheric 


flow over the sea surface. The flowfield cannot be strictly classified as a boundary laver 


= 
сл 


flow since non-uniform boundary conditions under constant fluctuation and changes 
almost alwavs exist [Ref. 32]. 

The lowest ten percent of the ABL is referred to as the surface laver and 1s 
tvpically 50 meters in height but may vary from 5-200 meters. The surface laver is 
characterized bv distinct variations in wind speed, nearly vertical heat and mass fluxes, 
and other meteorological fluctuations with height [Ref. 38: p. 256]. Fluctuations in wind 
variations, especially at sea may reach fifty percent of the mean wind speed 
[Ref. 31: p. 557]. The surface layer is also characterized by intense small scale turbu- 
lence which is mechanically produced by surface roughness or friction from waves on the 
ocean surface.6 The lowest portion of the surface layer is termed the roughness layer and 
is where the surface has the greatest effect on generating turbulence and influencing fluid 
motion. Mean flow in the roughness layer in largelv three dimensional and in- 
homogeneous in nature. Small scalc turbulence dominates flow near the surface. As al- 
titude increases, small scale turbulence decreases rapidly as large scale turbulence formed 
bv convective conditions in the atmosphere prevail. Several roughness heights above the 
surface, the mean flowfield becomes nearly horizontally homogeneous and two dimen- 
sional. Over sea, the mean surface laver and especially the roughness laver are more 
inhomogeneous and three dimensional than over land. 

2. Turbulence 

The roughness length zy is a measure of surface roughness as a function of the 
mean wind velocity and height above the surface. Surface laver turbulence length scale 
and intensity can be empirically determined from the aforementioned parameters. For 
typical open ocean roughness lengths in the range 10° = z, = 10 meters, turbulence 
intensities are on the order of 13 to 17 percent at a 10 meter elevation and mean wind 
speed of approximately 25 m sec [Ref. 39]. These high turbulence intensities are for high 
crosswind velocities and rough seas; extreme conditions which may be encountered. 
Tieleman states an average turbulence intensity of 6 to 12 percent to exist below a 20 
meter height above the ocean surface [Ref. 32]. 

From previous discussions, the turbulence length scale to body scale ratio is an 
important parameter for flowfield behavior. For the conditions stated in the preceding 
paragraph, the longitudinal turbulence length scales L, are in the range 80 « L, « 90 
meters (or 262 < L, < 295 feet). For a typical missile with a 1.1 foot diameter, the tur- 


bulence to missile diameter length scale ratio is about L,:L, ~ 280:1 and would have 


6 This small scale turbulence 1s larger than the missile length. 
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little effect on boundary laver development. Yet the possibility exists where boundary 
laver and vortex scale turbulence may be present due to the cascade effect from large 
scale turbulence and from crosswind interaction with the ship's superstructure, lending 
themselves to vortex and boundarv laver interactions. 

3. Crosswind 

A tvpical VLSAM launching with a 10g acceleration reaches an altitude of ap- 
proximately 56 feet travelling at a vertical velocity of V,~164 ft;sec. The VLSAM is still 
well within the surface layer environment (50 m) and subject to crosswind and turbu- 
lence effects. The preceding paragraph indicates turbulence length scales are relatively 
long in the surface laver. In such a case, the large scale flowfield or crosswind may dic- 
tate aerodvnamic effects on the VLSAM. As an example, if a ship travels at 20 knots 
and a mean wind speed of 20 m sec are combined with the VLSAM launch velocity, the 
resultant vector 15 at 191 ft’sec at an angle of 31° from the missile direction. The result 
is that the VLSAM flies at an effective angle of attack of 31? at 191 ft/sec, placing it in 
Regime III, the asymmetric vortex region, as it travels through the surface laver. 
[ ef. 3) 

A missile at a lower altitude, closer to the ship, will have a lower velocity and 
therefore it will be flying at an even greater effective angle of attack. Gregoriou states 
that a VLSAM that pitches over towards a target mav experience an effective angle of 
attack up to «x50? [Ref. 2]. These examples should illustrate the distinct possibility of 
asvmmetric vortex induced side forces on a VLSAM during launch phase. 

4. Ship Airwake. 

Changing flowfields and turbulence may be dictated by the ship’s hull and 
superstructure. This flowfield or ship airwake may increase the crosswind velocity, cause 
significant crosswind gradients, and non-uniformities in the flowfield increasing turbu- 
lence intensities while decreasing turbulence length scales. These effects occur in and 
around the vicinity of the ship which acts as a turbulence generator. Results of earlier 
discussions produced the following parameters in the marine surface laver: 


Mean crosswind velocity ОЕ ЗО (im Sec) 


tpemiencemtensie - 13< [p< 17 (percent) 


Turbulence length scale $0 « L,« 90 (meters) 
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These values mav change due to ship airwake effects in favor of higher turbulence in- 
tensities and smaller length scales, conditions prone to affect the missile boundary layer 
causing delaved separation. [Ref. 39] 

5. Additional launch considerations 

Other considerations which may influence the VLSAM at launch are dependent 
on the launch platform, the mechanics of a canister launch, and the aerodynamic effects 
of the the mussile controlling itself. These conditions are essential to the overall aero- 
dynanuc effects on the VLSAM as it exits the canister. 

The launch platform undergoes constant roll, pitch and yaw motion and trans- 
mits this motion to the VLSAM as it travels through the the guide mechanism in the 
canister. Increased severity of the ship’s motion adds complexity to the initial velocity 
vector as the VLSAM leaves the guide. 

The launch feature unique to a VLSAM is that the canister used for storage is 
also the launcher. Unlike a rail launched missile where motor exhaust dissipates into the 
surroundings, a canister or tube constrains a portion of the rocket motor gases resulting 
in a muzzle blast effect. As the VLSAM rocket motor ignites, most of the exhaust gases 
are typically vented via a plenum and uptake arrangement. Still, there 1s a significant 
amount of backflow which flows around the VLSAM as it leaves the canister. The 
rocket motor exhaust from the uptake may also interact with the flowfields about the 
nussile. The environment of a rocket motor ignition is severe in itself, with high mag- 
nitudes of heat, shocks, and vibration. [Ref. 40] 

The VLSAM itself may impose some form of flight control. aerodynamic or 
thrust vector, to maintain flight attitude during the launch phase. Should the nussile 
change flight attitude, the flowfield around it will be altered. 

These considerations are introduced to bring out the numerous factors which 
mav influence the VLSAM flowfield and cause in-plane and out-of-plane forces and 
moments during the launch phase. They will add complexity to a fully comprehensive 


flowheld analysis of the VLSAM which is bevond the scope of the present investigation. 


ПП. EXPERIMENT AND PROCEDURES 


A. APPARATUS 

The experimental investigation was conducted at the Naval Postgraduate School 
wind tunnel test facility. The VLSAM model was mounted on a six-component strain 
gage balance and sting support assembly. Turbulence generating grids and the test fa- 
cility steady-state data acquisition hardware/software complete the experimental appa- 
ratus discussed in this section. 

1. Wind Tunnel 

The horizontal-flow low-speed tunnel used for this experiment is located in the 
Naval Postgraduate School wind tunnel test facility and shown in Figure 5. Тһе wind 
tunnel is powered bv a one-hundred horsepower electric motor driving a three blade 
variable-pitch propeller via a four-speed Dodge truck transmission. Immediately fol- 
lowing the propeller blades are a set of stator blades to help straighten the flow. Turning 
Vanes installed at each corner, and two fine wire mesh screens placed at the entrance to 
the settling chamber smoothen the airflow and reduce turbulence. The settling chamber 
to test section contraction ratio 1s approximatelv 1071. A heavy wire mesh Screen 15 1n- 
stalled at the end of the test section diffuser to prevent foreign object damage to the fan 
section. [Ref. 41: pp.3-1 to 3-2] 

Тһе test section operates at atmospheric pressure and measures 45 inches by 32 
inches. Corner fillets house the test section lighting and help reduce svnergetic boundary 
laver effects at the wall interactions.? To counteract a reduction in freestream pressure 
due to boundarv laver growth. the walls of the test section are slightly diverged. Im- 
mediately downstream of the test section are breather stots around the perimeter of the 
wind tunnel. These slots allow air to enter and leave the tunnel unhindered. replenishing 
air lost through leaks and ensuring a uniform test section pressure. A reflection plane 
located in the center of the test section provides a horizontal surface to the flowhield and 
reduces the test section available height to 28 inches. A flush turntable at the center al- 
lows model the operator to change model pitch or angle of attack during wind tunnel 


operauon.[Ref. 41: pp. 3-3 to 3-4] 


7 The fillets reduce the area of the test section from 10 ft.? to 9.88 ft.2 Fillets are employed 
throughout the wind tunnel at wall interactions. 
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Figure 5. Naval Postgraduate School Low Speed Wind Tunnel: [Ref. 41: p. 3-7] 


Flow measurement components utilized with the wind tunnel are a dial 
thermometer, manometer, and pitot static tube. Temperature in the wind tunnel is 
measured bv a dial thermometer extending into the settling chamber. Dynamic pressure 
in the wind tunnel is measured by the static pressure difference between the test section 
and the settling chamber using a water filled manometer. Static pressure in the settling 
chamber is measured bv four static pressure taps, one on each wall. The test section has 
another four static pressure taps, one on each wall, upstream of the test section to pre- 
clude interference from the model. The pressure taps at each section are connected via 
a common manifold prior to feeding into the manometer. The manometer measures 
pressure differences in centimeters of water, which is the test section dynamic pressure, 


and is converted to the actual wind tunnel velocity by equation [5]. 


20 


U E 2 ° Πο e CM 11,0 ς 
B V 0.93 +p (5) 


ft 
sec 


Where: 
(и = measured velocity ( 
2.046 = conversion factor 
cm H,O = manometer reading in cm of H,O 
0.93 = calibration factor.8 


p = air density (πε 


A pressure transducer circuit, connected to the manometer displays test section dynamic 
pressure on a digital readout. Once calibrated, the digital manometer has a greater ac- 
curacy than the manometer and was used to measure the dynamic pressure. The pitot 
Static tube 1s read from. a aircraft airspeed indicator and is used only for test section ve- 
locity approximations when starting and stopping the wind tunnel. 
2. VSLAM Model 

The VLSAM model was designed to be representative of a cruciform tail-control 
missile with very low aspect ratio wings or dorsal fins. The model was fabricated from 
606] and 2024 aluminum alloy by Naval Postgraduate School personnel. The model ac- 
commodates the one-inch precision balance and sleeve used in the experiment and has 
internal reserve space for fitting pressure measurement instrumentation. Four sections 
form the major parts of the model. The body section is a hollow cylinder with locating 
pin attachment points for the balance. sleeve, wings and tails. The body mounts on the 
balance. or by an adapter for use without the balance to mount the model on the sup- 
port or sting. Body roll angle, with or without the balance, can be varied in 90° incre- 
ments. The nose piece attaches to the body forward of the balance. Nose roll angle may 
be vanied in 45° increments independent of the body. The nose 15 а pointed ogive shape 
with a short constant diameter section aft of the ogive which provides an interference 
free interface between the ogive shape and the body section. Four very low aspect ratio 
Wings with strakes comprise the cruciform wing section and are mounted equilaterally 
in fixed position along the model axis. Likewise, in the tail section, are four tail control 
fins which are fixed in postion at 0° incidence relative to the model axis mounted aft of 


each wing. All parts are rigidly connected to the model body by countersunk machine 


8 Total pressure correction for the settling chamber. 
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screws. The surface of the VLSAM model is free from proturberances and 1s bufled to 
a polished finish. Detailed dimensions of the VLSAM model are shown in Figure 6. 


A short summary of kev dimensions are listed below: 


Total lengeti = риа 
Base diameter = 1.75 in. 
Length diameter ratio = 13.06 
Ogive nose length = 4.0 in. 
Ogive diameter ratio = 2.29 
Wing span‘root chord = 3.13 1n./13.55 in. 
Tal span root chord — 5.50 1n./1.70 in. 


Center of pressure — 13.5 inches aft of nose tip (approx.) 
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Figure 6. Drawing of VLSAM model. 
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3. Balance 

Force and moment measurements were taken by a six-component, one-inch, 
Strain-gage precision balance. The balance was on loan from the NASA-Ames Research 
Center under the Navv-NASA Joint Institute of Aeronautics. Balance calibration was 
performed bv NASA-Ames. Calibration data and data conversion values were also 
provided and are listed in Appendix A [Ref. 42]. A wheatstone bridge circuit capable 
of measuring positive and negative strain comprises each channel and is oriented to 
measure strain in a chosen direction. The directions relative to the balance axis are 
normal, side, and axial with one gage measuring rolling moment. Channel capacities 
and accuracies are listed in Table 1. Individual channel accuracies are based on full 


scale (maximum) load. 


Table I. BALANCE CHANNELS: [Ref 42]. 


| Normal Force | Force BEIM lbs. | | 00856 . | | 00856 . 
= [езе [вы 
Ls [omae | oem L| oem 


A machined sleeve provides a close tolerance fit between the balance and the 










interior of the VLSAM model. The model is seated to the balance by locating pins such 
that the balance center coincides with the approximate center of pressure of the model. 
Since the balance supports the model, it provides the mechanical interface between the 
model and sting support. 
4. Mfodel/Balance Support 

The VLSAM model and balance support, also called the sting mount, 1s rigidly 
fixed in the wind tunnel test section bv the reflection plane turntable at the base and an 
aluminum reinforced clear plexiglas section at the top. The sting mount changes model 


pitch in the horizontal direction with the pivot point located to coincide with the ap- 
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proximate centers of the VLSAM model and balance. Heavy duty construction in 
stainless steel and aluminum keeps deflections and vibrations to a minimum. The pre- 
sented cross-sectional area of the sting mount to the test section flow is minimal for all 
pitch angles to reduce blockage effects. The reflection plane turntable rests on a 
heavy-weight pedestal to prevent deflections and vibration. The entire assembly rotates 
via a chain gear drive powered by an electric motor. Pitch and pitch rate are set by the 
operator on the motor controller, or by a relay actuator controlled by the data acquisi- 


tion system. A ring scale on the pedestal indicates model pitch in 1/10 degree increments. 





Figure 7. Photograph of VLSAM model in the test section. 
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Figure 8. Planview of VLSAM model in the test section. 


5. Turbulence Grids 

The turbulence grids were designed for use in the low speed wind tunnel to 
generate turbulence in Varving intensities and length scales. Each of the four grids or 
screens 1s mounted in a wooden frame and placed 73 inches forward of the model sting 
mount pivot point. Tbe specifications for the four grids are listed in Table 2. Three of 
the four grids are square-mesh square-bar and fabricated from wood. They employ bi- 
plane construction which generates a nearly isotropic homogeneous turbulence. The 
fourth grid is also square mesh, but was fabricated from round wire. Two wire screens 
are attached together for a fine mesh. Turbulence intensities and length scale for each 
of the grids are shown in Figure 9 and Figure 10. These grids were constructed by 


Roane who also measured the turbulence intensities and length scales. 


Table 2. СКІ) ВЕСИ СИМО ОИ) 


















Grid Nlesh Width (in.) Material 


Bar Diameter (in.) | Mesh/Diameter 














22 

20 

18 - D GRID 1 
^ GRID 2 

16 9 GRID 3 
* GRID 4 


— 
+ 


о 


TURBULENCE INTENSITY (PERCENT) 


80 90 


0 10 | 20 30 40 50 60 70 
DISTANCE FROM GRID (INCHES) 


Figure 9. Grid Turbulence Intensity: The dashed line indicates the model pivot 


axis. [Ref. 3: pp. 45] 
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Figure 10. Grid Turbulence Length Scales: The dashed line indicates the model 


pivot axis. [Ref. 3: pp. 48] 


Table 3. TURBULENCE PARAMETERS: At the model pivot axis. 
[Ref. 3: pp. 44-49] 


Grid Intensity Length Scale Turbulence/ Dynamic Pressure 
(percent) (in.) Model Dia. (Ib/ft?) 
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The dynamic pressures shown in Table 3 are the true dvnamic pressure in the test sec- 
tion from Roane [Ref. 3: pp. 44-49].. The grid generated turbulence length scales to 
missile length ratios at the test section are smaller than 1:12 which suggests that the 
turbulence effects on the flowfield and vortices may be due to boundary laver scale 
and. or vortex scale turbulence. The turbulence length scales to body diameter ratios are 
on the order of and smaller than one as seen in Table 3. Conclusions on the effects of 
the grid generated turbulence with regards to changing length scales at constant intensity 
or changing intensities with constant length scale are not possible with the present grid 


turbulence parameters. 
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Figure 11. Turbulence Grids 2, 3, and 4: Clockwise from lower left hand corner; 
Grid 4, Grid 3, Grid 2. 


28 


| 
| 
| 
| 
| 
J 









® w 
«ἲ за : x hs Ced v. ' 
> iex š AN * es LS SESS ОНА 5 
Ж M RR š 5% ? wn i 
v UNES $i се s 
инна ш 
Е 3 7! 
4 - X 3 
e 
x i 3 3 
Ν΄ ο E i 
K asa 4 Жы ç * 


1 










ехе 
> 





] 


© WRATH 





Г] tet PR ee 
хе s 
ñ. | 
& 





ο red nd Cán nnd Lad uml 


z ^ 
$ 
%” 2, 
" Р 
- Ë s =s xd P 22e в š κκ... 4 
4 = " s 4 
ЖУ 
> z 
> АС 
š > | 
EN L^ 
2 ‘fe 
^ XY «e. ^, 9 n ” с РУ Е d % 
X T 2 et 2 /ж S Y * ὦ C 
dd 
| E. 
^ 
, 
ü 5 
° * ως. t 
© 
Р 





Figure 12. Turbulence Grid 1l installed in the wind tunnel. 


6. Data Acquisition Hardware 

The data acquisition hardware consists of test facility components to operate 
the precision balance in the wind tunnel. Each strain-gage in the balance is supported 
by a signal conditioning circuit which supplies excitation voltage and allows zeroing and 
calibration of the balance circuit. Each balance channel from the signal conditioner unit 
is fed into a Hewlett-Packard relay multiplexer which sequences through each channel 
as required by the data acquisition software. Prior to being measured, the balance signal 
is fed through a 1000 gain low-noise amplifier to improved resolution for small magni- 
tude signals. A Hewlett-Packard digital multimeter converts the analog voltage signal 
into a digital signal for use by the data acquisition program. An IBM-AT microcom- 
puter executes the data acquisition software and stores the collected data. Other items 
utilized in the data acquisition process are the wind tunnel measurement apparatus dis- 


cussed earlier. Data acquisition hardware components are shown in Figure 13. 
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Figure 13. Data Acquisition Hardware 
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7. Data Acquisition Software 

The data acquisition software consists of a shell program to control the 
Hewlett-Packard instruments and a data conversion program to convert the raw balance 
readings into a usable form and store them. All readings are time averaged and represent 
steady state measurements. The initial balance reading in volts is manipulated with the 
balance calibration constants to give force values in pounds and moment values in 
foot-pounds. Another program translates the force and moment values into a coefficient 
form accounting for test conditions and corrections in the wind tunnel. The data ac- 


quisition programs are listed in Appendix C and Appendix D. 


B. EXPERIMENTAL CONDITIONS 
The introduction discussed manv variables which could affect flow separation and 
vortex structure at high angles of attack. To simplify the various outcomes of the test 
due to numerous perturbations, the following variables were kept as constant as exper- 
imental conditions would allow: 
]. Revnolds number 
2. Test section dynamic pressure. 
3. Test section velocity. 
4, Nose geometry and roll angle. 
The experimental reference pressures used in the turbulence mapping by Roane 
[Ref. 5] were duplicated in this experiment in order to use the turbulence grid length 
scales and intensities shown in Table 3, Reference dynamic pressures for the grid and 
no grid conditions are listed below. These reference pressures ensure a sinular test section 
Velocity since the total pressure difference measured by manometer 1s in error due to the 
disturbances caused by the screens. Since actual dynamic pressure in the test section 15 
determined bv the test section velocity, the velocities in the test section from Roane’s 
turbulence mapping were used to calculate the actual dynamic pressures. These refer- 
ence pressures are: 
No GRID: a, = 7.20 (ст. Н.О) 
GRID 1: дк = 10.00 (ст. Н,О) 
GRID 2: q, = 10.00 (ст. Н,О) 
GRID 3: qg = 10.00 (cm. H,O) 
GRID 4: αι = 5.40 (ст. Н,О) 
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The dvnamic pressure conditions for the experiment would give a subcritical Reynolds 
number of Re,xl.1X10* and actual dynamic pressures in the test section shown in 
Table 3. Test condiuionsdorandpsdusiunstaredisted Ino ppo MH ME NOCHES ου 
were altered to make comparisons: 

]. Angle of attack. 

2. Wing tail configuration. 

3. Afterbodv roll angle. 

Three diflerent VLSAM model afterbodv configurations were chosen to compare the 

effects of turbulence on the VLSAM model with and without wings and at two different 
roll angles. 


1. BODY A - Wings and tails at zero roll angle in a “+” configuration. 


гә 


BODY B - No wings or tails. Body roll angle set at $4 — 45? . 


3. BODY C - Wings and tails at óg — 45^ roll angle in a ^N" configuration. - 
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Figure 14. Body configurations and reference system. 


In all cases. as the afterbody roll angle changes, the nose roll angle remains fixed to the 
reference plane. BODY A and BODY C configurations were chosen to be certain that 
measurements would indicate the presence of vortex asymmetry by any side force mag- 
nitude. The wings and tails in the “+” and “X” roll angles are the only wing attitudes 
where, in the absence of asymmetric vortices, zero side forces, rolling moments, and 
vawing moments occur. [Ref. 43] 

The goal of the experiment was to measure the effects of the grid generated turbu- 
lence on the high angle of attack aerodynamics on the VLSAM model, in particular the 


affect on side force magnitude and direction. 


C. EXPERIMENTAL PROCEDURE 
1. Balance Calibration 

The precision balance was calibrated by the NASA Ames Calibration Labora- 
tory which determined balance constants for the data acquisition program and gave in- 
dividual channel or circuit accuracy. The balance was fixed to the model'balance support 
in the test section without the sleeve and VLSAM model. After all electrical connections 
were made and checked the data acquisition hardware was energized and balance 
excitation voltage for each channel was set to 5.00 + 0.003 DC volts via the signal con- 
ditioner. The entire data acquisition svstem was allowed to warm up and settle out for 
one hour. Next, the amplifier was zeroed at a 1000 gain setting followed by zeroing each 
balance circuit to 0 + 200 nanovolts across the wheatstone bridge. The balance sleeve 
was fitted along with a cradle for static weights. The data acquisition software was run 
to obtain the initial tare readings. Finallv calibrated weights were added to the cradle and 
measured by the balance to compare readings with balance calibration specifications. 
Normal and side force readings were found to be accurate to within 0.05% of full scale 
readings. After this initial test was completed, the signal conditioner was never de- 
energized to maintain the excitation voltages to the balance circuits. 

2. Prerun Calibration and Test 

Prior to each run, the data acquisition system was fully energized and allowed 
to settle out for half an hour. During the settling period, the data acquisition system 
was cycled repeatedly by the software to bring the relay multiplexer, digital multimeter, 
and amplifier to operating temperature conditions. Next, the static weight test was re- 
peated with the mounted VLSAM model, to ensure the balance circuits were remaining 
within calibration accuracy. If not, the model was removed, the balance circuits were 


zeroed. and the static weight test procedure completed without and with the model. 
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3. Data Collection 
When the prerun calibration results were satisfactory, the wind tunnel was started and 
test section reference pressure set at 110% of the data run reference pressure. The model 
pitch or angle of attack was slowly cycled from -10° to 110° and back at approximately 
2° per second. This procedure tested the safety of the support assembly with the balance 
and the model and let the mechanical interfaces settle out. The pitch angles where 
noticable model vibration occurred were noted. 

The wind tunnel was put into a no wind condition (zero dynamic pressure) with 
the motor running and the model set at -5° pitch.9 The data acquisition software re- 
corded the initial tare readings. The intial settling chamber temperature and barometric 
pressure were recorded. The wind tunnel was brought up to operating reference pressure 
and the actual data acquisition began. Missile angle of attack was manually set and 
entered into the data acquisition program for each set of readings. Test section reference 
pressure was kept constant as the model angle of attack increased adding blockage. The 
microcomputer reduced the balance readings to force and moment values and stored the 
data on disk. At the end of each run, reference pressure was set at zero and several 
readings recorded to measure anv shift between initial and final readings. 

The force and moment data are translated into coefficient form by another 
program which also corrects for blockage and test conditions. Data were graphed for 
comparative analvsis and graphs are found in the results section. 

4. Preliminary Runs 

Previous experiments showed that variations in nose roll angle would alter the 
asymmetric vortex structure, thus changing side force direction and magnitude. Non- 
repeatabilitv of data points is another area often encountered in earlier invesugations. 
Prior to conducting the actual experimental data runs. preliminary runs were made to 
determine which nose roll angle produced the largest side force magnitude. Data runs 
were also conducted to examine the repeatability of data without changing test condi- 
tions. The preliminary runs would also indicate data trends and would give some degree 
of confidence in the data collected in the subsequent runs. 

Preliminary test runs were short and durations ran approximately 30 minutes. 
Temperature increases were small, from 2^F to 5^F, giving a smaller temperature gradi- 
ent and allowing test section air density to remain in a narrow range. When the exper- 


iment proceeded to the actual test runs, test time durations tvpicallv lasted 4.5 hours. 


9 This conditon will be referred to as the model's zero condition. 
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In order to keep test section air density as uniform as possible, the test procedure was 
interrupted at least once during the run cycle to cool the air in the wind tunnel. Settling 
chamber temperatures rose about 10°F to 20°F with larger increases for coarser turbu- 
lence generating grids. When the the settling chamber temperature increased to 20°F 
above the start temperature, the wind tunnel was cooled down by air exchange. Once 
the tunnel had cooled down to within 5°F of the start temperature, the run continued 
at 2° below the interrupt angle of attack to overlap readings. 

Other test conditions that were beyond the control of the experiment include the 
turbulence and vibrations produced by the wind tunnel itself during operation and vi- 


brations and deflections of the model balance support. 


D. EXPERIMENTAL CORRECTIONS 
Model blockage corrections for changing angle of attack were made by equations 
(6], [7]. and [8], from [Ref. 44] 


а = quill + 22) {6} 
U z UA(1 4 £) Ш 
where: 
q = dynamic pressure (ib, ft?) 
cneasuseg referemee pressure (lIb/ft2) 
U = horizontal velocity (ft/sec) 
U,, = measured horizontal velocity (ft/sec) 
€ = blockage factor 
and: 


] model frontal area (8) 
E аа ecc C 
34 test section area 


£ = 
Model frontal area varies for each body configuration and angle of attack with minimum 
blockage at 0° and maximum blockage at 90° as shown in Figure 15. 
Blockage equations for each body configuration were developed as a function of 
model angle of attack. These equations are shown in Figure 15 and Table 4, and are 


implemented in the data conversion program listed in Appendix E. 
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Table 4. VLSAM MODEL BLOCKAGE 


Minimum Blockage Maximum Blockage 
(percent) (percent) 


Configuration 
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Figure 15. Blockage factors 


36 


HI. RESULTS 


Several important factors must be kept in mind as the results are presented. These 
factors are conclusions from previous high angle of attack studies on slender bodies in 
low turbulence environments and may apply to these results. The nose fineness ratio of 
the VLSAM model 1s 2.29. Results bv Keener [Ref. 11] dictate that the vortex cores from 
the nose at fineness ratio 3.5 should not be as crowded together as in higher fineness 
ratio noses at 5.0. Since the nose is pointed, no blunt nose effects are expected. From 
equation [2]. vortices should become asymmetric when the angle of attack is a,y=24° for 
the given apex angle of the VLSAM model. Flow separation throughout the length of 
the nose should remain laminar at all angles of attack for Re, x 3.0x10: [Ref. II]. Con- 
ditions for all data runs are stated in the previous chapter and listed in Appendix B are 
stated in the previous chapter. All data represent steady-state conditions and are ex- 


pressed in coefficient form. 


A. VARIATIONS IN NOSE ROLL ANGLE 

Previous research showed that variations in nose geometry and roll angle altered the 
Vortex structure shedding from the nose and afterbody at high angles of attack. The 
objective of the first set of preliminary runs was to determine the nose geometry which 
produced the highest side force magnitude. Since the model used in the experiment used 
onlv one nose, variations in nose roll angle became the single variable affecting nose 
geometry. The nose with the highest side force magnitude would be fixed in the same 
position regardless of afterbodv configuration for all subsequent runs. Both sets of runs 
were Conducted under identical conditions, Re, = 1.33x105 and test section. dynamic 
pressure at q = 22.2 lbs ft?. Angle of attack was varied from 0° to 90° in 10° and 5° 
increments. Model configuration was BODY A and no grids were used. Nose roll angle 
was changed in 45° increments for each run. Results for each set of runs are plotted in 
Figure 16 and Figure 17. 

In both cases, changes in nose roll angle alone produced variations in side force 
magnitude and direction. The plots also show a some repeatability of the data points. 
Nose eight, in both run sets, produced the largest side force magnitude in the negative 
direction. The angle of attack when side force magnitude became noticeable at oz: 30* 
was the same for all nose roll angles within each set of runs, which 15 also true for the 


angle of attack «x $0? where side force magnitude decreases to negligible values. 
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Figure 16. Side force variations with nose roll angle: Runs W04A101 to WOASOI. 
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These side force measurements suggest that vortices become asymmetric at ο, 20° апі 
end сз which is true for all roll angles. The angle of attack increments are large. 
Smaller increments would give better observations of data trends and all subsequent data 


collection runs are in one degree increments. 


B. DAIA REPEAIABILITY 


A recurring problem dealing with side force measurements from asynumetric vortices 
15 the non-repeatability of data except in Keener (Ref. 23] where repeatability of results 
was good. In this experiment, all readings are an average of ten separate measurements 
taken within a 40 second interval, so any signs of repeatability could indicate that the 
data 1s statistically accurate. As shown in Figure 16 and Figure 17, the nose roll angle 
which gave the maximum side force magnitude for all roll angles 1s nose eight, which 1s 
true for both run sets. The objective of the repeatability runs was to determine if any 
repeatability existed, and if not, what degree of fluctuation existed. After prerun cah- 
brations, the model was pitched to 20? angle of attack from 0° at test dynamic pressure. 
Readings were taken at 20° then the dynamic pressure reduced to zero and readings 


taken again. The wind was turned on and off for readings as the model remained at a 
fixed angle of attack. 
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Figure 18. Data Repeatability Runs 
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After several readings. the model was slewed to the the zero condition for a final cali- 
bration check. This procedure was repeated for 45°. 70°, and 90° angles of attack and 
all results are plotted in Figure 18. Variation in side force magnitude 1s obvious at 
20°, 45°, and 70°. Measurements at 90° show negligible fluctuation. This agrees with 
accepted side force fluctuation phenomenon. Staustically, the mean value of anv data, 
will be improved if an increasing population of readings are taken. Time and record 
keeping limitations allowed only ten readings per channel to be recorded. Some degree 
of experimental accuracy mav be lost with the data fluctuation, but the data point trends 
are useful in making comparisons between different run conditions. Comparing 
Figure 16 and Figure 17, the side force follows a general trend which will be useful for 


comparing the data runs discussed in the following. sections. 


C. BODY IN ISOLATION 
The baseline test run was made for BODY B with no grids shown 1n Figure 19. 
The asvmmetric vortices begin at &4,526" which concurs with the value o,:~24" from 


equation [2]. The unsteady vortices begin to form at o,,280° . Between these two ап- 
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Fieure 19. КОВ802: Side Force Coefficient. 
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cles of attack the side force changes magnitude but not direction. The data points 
scatter suggests a strong possibility of fluctuations in the steady state side force readings. 
The data trend shows the maximum side force occurring at ax56? and C, 2 4.1 then ta- 
pers off to nearly C, — O after which it increases back to C, = —2.0 then tapers off into 
the unsteadv vortex region. Since there are no wings or tails, the side force is due en- 
tirelv from the nose and afterbody vortices. The smaller angle of attack interval in Fig- 
ure 19 does not show any data smoothness when compared to Figure 16 and 
Figure 17, which were run at larger increments. The data trends are similar which illus- 
trates the dominance of the nose generated vortices. The side force magnitudes in 
Figure 16 and Figure 17, are just over 30% of those in Figure 19 which may be due to 
wing vortex effects discussed later in the results section.10 

The BODY B configuration with Grid 1 generated turbulence is shown in 
Figure 20. The onset of asymmetric vortices is delayed to o,y=42° and unsteady asym- 
metric vortices a-.=82° . In general the side force direction remains the same but there 
is a considerable reduction in magnitude by about 77%. One possible cause of this de- 
crease in side force 1s the change in positioning of the asymmetric vortices caused by the 
turbulence. The turbulence length scale to missile diameter ratio 1s 1.08 and can be 
classified as vortex. scale. turbulence which may suppress the asymmetry of the nose 
vortex pair. Perhaps this 15 the reason behind the delayed angle of attack for the onset 
of asvmmetric vortices. Another reason may be due to turbulence intensity since 
GRID 1 turbulence intensity is the highest. 

Comparing the results from the remaining grids in Figure 21 thru Figure 23, the 
trends are similar to GRID I and the vortex asymmetry seems to be depressed by longer 
turbulence scales. For GRID 4 where the turbulence length scale 1s the smallest and 
side force magnitude is the highest, xay occurs at an earlier angle of attack. The data 
trend in Figure 23 has a closer similarity to Figure 19, although side force magnitude 
is reduced by 55%. Turbulence intensity does not seem to factor as much as the length 
scale to control asymmetric vortices since turbulence of any intensity might be expected 
to eliminate asymmetric vortex effects and especially in the case of GRID 1 with the 


highest turbulence intensities, yet GRID 1 results show the side force is not eliminated. 


10 The axisymmetric wing and tail arrangement of BODY A and BODY B should not con- 
tribute any side force lift from the lifting surfaces. 
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Figure 20. КІВ801: Side Force Coefficient. 
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Figure 21.  R2B80E: Side Force Coefficient. 
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Figure 23. R4B801: Side Force Coefficient. 
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A correlation between turbulence length scale and turbulence intensities to the vortices 
is not possible since the turbulence intensities for the grids are not equal. 

In previous research, Wardlaw produced 2% turbulence intensity by a wire mesh 
grid [Ref. 12].11 His results showed that side force distribution was not effected by the 
turbulence and side force magnitude increased by 10%. If the turbulence length scales 
were small enough to affect the boundary layer, the side force magnitude should de- 
crease. On the other hand, larger turbulence length scales should also suppress side 
forces as indicated by the results in this experiment. Since the small scale low intensity 
turbulence from GRID: 4 does not appear to suppress side force magnitude, results from 
increasing GRID 4 turbulence intensitv could settle whether turbulence intensitv con- 
tributes anvthing significant to the vortex structures. 

In another investigation bx Clark, he showed that a single rod place upstream from 
the stagnation point on an ogive nose may have the same effect as grid generated tur- 
bulence. He injected a thin stream of bubbles into the flowfield upstream of the model 
which forced flow separation into steadiness [Ref. 45].12 His results indicated that when 
bubbles were mjected at the pomt oft the nose, the vortices were forced to гета“ 
fixed separation condition. When Clark injected the bubble stream on the afterbodv, 
there was no effect on the vortex structure or side force magnitude. He concluded that 
the wake from the bubble probe, rather than the bubbles themselves, forced the flow 
separation into steadiness. 

Conjecture at this point would indicate that turbulence length scales of vortex size 
mav have a greater influence on the vortex structure than turbulence иеле ee 
smaller scale low intensity turbulence does not seem to attenuate the asymmetric vortex 
Structure since maximum side force magnitudes are nearly equal, but does seem to cause 


a steadier flow evidenced bv the data scatter. 


D. WINGS AND STRAKES WITHOUT TURBULENCE. 
The addition of wings and strakes in the BODY A and BODY C configurations 
show changing trends in side force direction and magnitude.13 BODY A in Figure 24 


shows larger side force, a lower angle of attack for the onset of asymmetric vortices, 


11 Ве, = 1.55 ООСН που. 
12 Ке, = 0.4х10°. No turbulence length scale or intensity given. 


13 The following discussion centers on the wings and strakes since the tails are not expected 
to contribute any significant effects. 
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larger side force magnitudes by 1290 when compared to BODY В іп Figure 19, and 
greater flow unsteadiness seen in the data scatter. BODY C in Figure 25, while showing 
an equal degree of unsteadiness, has even lower side force magnitudes than BODY B 
and pronounced changes in the direction of the side forces seen from the positive side 
force values between «%18° and «=50°. These results suggest that the addition of wings 
do not significantly alter side force magnitude. The nose vortices continue to dominate, 
but the wings and strakes do influence side force direction which could be caused by 
other flow mechanisms. 

The afterbody vortices of the BODY B configuration are expected to contribute side 
force to the overall magnitude which raises questions as to what, if anv, are the position 
of afterbodv vortices in the BODY A and BODY C configurations. The symmetry of 
BODY A and BODY C configurations are not expected to contribute any side force lift 
from the wings and tail in the absence of asymmetric vortices. Since side force magni- 
tude and direction trends are generally similar but not alike for all cases, there must be 
other mechanisms involved froin the vortices generated by the leading edges of the 
strakes and wings. 

The first possible effect is the strength and positioning of the wing and strake 
vortices. The nose vortices are alwavs expected to dominate any downstream vortices 
and may drive the wing and strake vortices to asymmetry contributing to the overall side 
force magnitude. These wing and strake vortices could be stronger at inducing side 
forces than the afterbody vortices and may be the reason why BODY A has larger side 
force magnitudes than BODY B. The second effect may be due to the strake and wing 
vortices pushing the nose vortices further awav from the body causing a weaker asvm- 
metric Vortex effect and reducing the overall side force magnitude which may explain 
why BODY C has smaller side force magnitudes than BODY B without turbulence. It 
mav also be a combination of the two and could be dependent on the angle of attack. 
The third effect may be due to the nose and afterbodv vortices causing side force lift on 
the leeward wing. As the nose vortex pair becomes asymmetric, the larger vortex could 
roll up and over the leeward wing causing a side force from vortex lift, a phenomenon 
enhanced bv the strake and low aspect ratio wing. Side force vortex lift could also 
contribute to the rolling moment. Finallv, the wings and strakes, depending on roll an- 
gle. mav act as splitter plates which form a barrier between the vortex pairs preventing 


or attenuating an interaction. 


High angle of attack experiments by Deane [Ref. 27] and Jorgensen [Ref. 46] with 
monoplane configurations show no increase in side force magnitudes when compared to 
the body in isolation. Further studies with monoplane wing configurations at 0° and 
90° roll angles would be invaluable to study these possible effects. The disposition and 
Strength of the asymmetric vortices would have to be quantified to determine how the 


actual vortex interactions occur. 


E. WINGS AND STRAKES WITH TURBULENCE. 

The effects of adding turbulence to BODY A and BODY C were similar. In both 
cases Figure 26 through Figure 33, the longer turbulence scales suppressed the side 
force magnitudes, but. for BODY C the data show more unsteadiness with the smaller 
scale turbulence. as opposed to the BODY B results. The larger scale turbulence delavs 
the angle of attack for the onset of asvmmetric vortices for BODY C to o,,x23? The 
results for BODY A are again close to BODY B with the onset of asymmetric vortices 
Occumane earlier at x4-22^ for the no turbulence condition. The effects of the larger 
scale turbulence in delaving z4, are apparent when GRID 1 results in Figure 26 and 
Figure 30 are compared to GRID 4 results in Figure 29 and Figure 33. Changes in side 
force magnitude are mixed although the general trend shows the lower side force mag- 
nitude at longer turbulence length scales except for BODY A and GRID 3 m 
Figure 28 which has a higher side force magnitude. The side force directional tendencies 
are also suppressed suggesting longer turbulence length scales brings more steadiness to 
the asymmetric vortices. Turbulence intensity effects from the graphs show that the 
higher intensities suppress side force magnitude but do not seem to change flow sepa- 
ration unsteadiness because of the data point fluctuation. The general trend of the data 
points are sinular for all cases examined. especially when comparing each separate body 
configuration to each other. These trends indicate the dominating influence of the nose 
vortices. The side force differences between the body configurations could be attested to 
the presence of the wings and their roll angle. 

Vortex scale turbulence may be responsible for changing the position of the vortices 
which explains the changing direction and magnitude of the side forces which 1s true for 
the case of BODY B with GRID 1. Here the vortex scale turbulence may contribute to 
the general unsteadiness of the vortex structure evidenced by the scattered data points. 
This contribution to unsteadiness can be compared to a missile experiencing unstead- 


iness with comparable turbulence length scales. Boundary laver scale turbulence mav 
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exist due to the cascade effect but does not seem to factor into the experimental results 


since side force magnitudes were always larger for smaller scale turbulence. 


F. NORMAL FORCE RESULTS 

The normal force results are plotted in Figures 33 through 47 and show the charac- 
teristic trend of a steady increase to a maximum where the values level off. The normal 
force is largest for the BODY A configuration and the smallest for BODY B. The 
maximum normal force value of the BODY А configuration is C, — 32.3 which 15 22% 
higher than BODY C at C, = 26.4 and 101% higher than BODY B at C, = 16.1 The 
normal force curve for BODY A is much steeper and levels off at the higher value and 
lower angle of attack than the other configurations. This may be due to the wings 
producing more lift at lower angles of attack and more drag at higher angles of attack 
than for BODY C. The addition of grid generated turbulence causes a decrease in slope 
of the normal force curves. The curves without turbulence show a sharp turn at a max- 
imum value where the curve levels off. The curves become rounded with the turbulence. 
There does not appear to be any significant effect on the overall magnitude of the 
normal force in the presence of turbulence. The variation in turbulence length scale and 
intensity do not appear to affect the normal force values as in the results for side forces. 
There is no noticeable fluctuations in the trend of the data points since the normal force 


does not fluctuate with or without the presence of turbulence. 


G. SIDE FORCE TO NORMAL FORCE RATIO 

The side force to normal force ratio 1s important to slender bodies with monoplane 
wings such as cruise missiles and aircraft. As this ratio increases, the slender body must 
be able to cope with the increasing side forces, thus requiring large amounts of vaw 
control. Since yaw control does not equal lift control due to the available control surface 
area in monoplane configurations, side forces can be a problem. Cruciform missiles 
should be able to cope with high side force to normal force ratios since their wing and 
tail arrangement gives equal lift capability in the pitch and yaw directions. 

Since the normal force components are not as effected by the addition of grid gen- 
erated turbulence, the side force to normal force ratios should follow the side force 
trends. The lowest normal force values of all the tested body configurations was 
BODY B which also had a substantial side force value, and therefore should be expected 
to have the highest side force to normal force ratio which is the case as shown in the 
plots. BODY B has a maximum ratio of C, / C, = —0.369 which is 112% higher than 
BODY A at C,/ C, =—0.174 and 367% higher than BODY Е зме г = —0.079 . 


po 


The side force to normal force ratios for the turbulence conditions in. Figures 50 
through 63. should follow the side force (with turbulence) trends since the normal force 
components are not affected to the same degree as the side force components. The 
C, ' C, plots show some resemblance to the trends in the side force plots. The wide 
degree of scatter. especially for the turbulence from GRID 4 is representative of the 
side force trends. The cruciform nussile should be able to cope with higher side force to 
normal force ratios due to it’s wing and tail configuration. From these results, the 
winged configurations, BODY A and BODY C have smaller side force ratios than the 
body in isolation of BODY B. Although the ratio is small, this may be misleading since 
the BODY A and BODY C configurations must be able to cope with higher normal and 


side forces than BODY B which can be easily accomplished by their cruciform wing 
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Figure 34. .ROBS02: Normal Force Coefficient. 
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Figure 39.  ROASQ3: Normal Force Coefficient. 
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Figure 42. R2A801: Normal Force Coefficient. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The VLSAM in the marine environment may encounter turbulence during the 
launch phase at high angles of attack. The turbulence intensities in the marine envi- 
ronment are higher than produced in the experiment, as are the turbulence length scales 
which are on on the order of LL, x 280:1 for a typical VLSAM. Cascading and or ship 
airwake mav produce smaller low intensity turbulence on the order of missile diameter 
or vortex size. These turbulence scales relative to the model diameter тау Бе produced 
іп the wind tunnel. 

Subsonic wind tunnel investigations were performed on a cruciform wing missile 
model with different grid generated turbulence length scales and intensities. The model 
was tested with three different body configurations with and without wings at 0° to 90° 
angle of attack. The results from the experimental runs with no turbulence yield the 
following conclusions: 

1. Side force magnitude and direction change with vanations in nose roll angle. 

2. There is a general trend in the repeatabihty of data points and to some degree for 
all data points (for the steadv state condition) which suggests the dominating influence 
of the nose generated vortices. 

3. Strake and wing vortices may play an important role to influence the disposition 
of nose generated vortices. While the low pressure on the top of the wing may pull the 
nose vortices closer to the body, leading edge vortices from strakes and wings may push 
them awav. 

J. The maximum side force-to-normal force ratio was observed for the bodv in 1so- 
lation configuration, although normal and side force values were actually lower than the 
winged configurations. 

The results of the turbulence experiments indicate that turbulence effects for side 
forces may vary with length scale and intensity. Two turbulence length scales, boundary 
laver scale turbulence and vortex scale turbulence, are particularly important as the 
former may have large effect on flow separation mechanisms which generate the asym- 
metric vortices and the latter on the disposition of the vortices in the flowfield. Turbu- 
lence results conclude: 

l. Longer turbulence length scales decrease side force magnitude and possibly any 


flow unsteadiness in the asymmetric vortex structure. 
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2. The angle of attack for the onset of asymmetric vortices is postponed bv longer 
turbulence length scales. Since these length scales are quite large when compared to the 
boundary layer, the asymmetric vortex structure may be suppressed by the interaction 
between the inviscid vortex flowfield and the vortex scale free stream turbulence. 

3. Higher turbulence intensities cause less flow unsteadiness for missiles with wings. 

4. Turbulence had little effect on the normal force values. Maximum lift values were 
nearly identical and there was no fluctuation in normal force values for all turbulence 
conditions. 

Further recommendations for a continued study of wing and turbulence interactions 
include: 

1. Perform flow visualization, hotwire, surface flow, laser doppler, and water tunnel 
tests to quantify the various flowfields about the missile model at high angles of attack. 

2. Test monoplane wing configurations at 0° and 90° roll. 

3. Compare the effects of turbulence length scales at the same intensities and with 
longer and shorter turbulence length scales than those used in this experiment. 

4. Compare the effects of turbulence intensities at the same length scale and with 
higher intensities than those used in this experiment. 

5. Model the rocket motor plume to study the interaction with flowfields about the 


missile. 
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NDEP'ENDIN-'A. 


BALANCE CALIBRATION CONSTANTS 


BALANCE CALIBRATION CA. DATE 7247 COMP JATE 7287 
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ЗА -АРА-ТТТ МАХ 10%0 ІНч5 X 5345Е СА. SHUNT CAL ROS 
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м2 400.00 600.00 359. 0.1657 100K «626 
A 100.00 100.20 175. S)K &618 
51 200.00 200-00 350. 0.1375 0 OK 4623 
S2 200.00 200.00 352. 0.1375 10 0K «597 
RM 21.00 20.83 175. SOK 4623 
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ЗА K_POSC1) Қ Р05С2) К МЕСС1) K_N=GC2) DEV АС. 
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М2 %.7211Е-02 -1.7015Е-08 %8.7763Е-02 8.9153Е-09 0.196 0.059 
A. 1.4309E-02 -7.1962€E-10 1.4290E-02 -1.3322E-09 0.115 0.115 
51 3.1309Е-02 -3.8153Е-08 3.20273Е-02 -8.9316Е-09 -0.263 0.13? 
$2 3.0366c-02 -3.8607TE-08 3.1167Е-02 -7.251ТЕ-09 0.3255 0.153 
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Fieure 64. 
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Balance Calibration Constants: 
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N1/N2- к-1.92257:-22 
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[Ref. 42] 
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APPENDIX B. RUN MATRIX 


The run matrix lists all runs conducted for this thesis research. All runs were con- 
ducted bv the author in the Naval Postgraduate School Wind Tunnel Test Facility. Not 
all of the runs were discussed in the results section and were essential trial runs to fa- 
miharize the author with the data acquisition system and observe the general trend of 
the data. 


These runs are listed bv a sequence code as in: 
ROBS02 
where: 
1. The first letter 1s a run code. 
"W* - Nose roll angle runs. 
"E” - Data repeatability run. 


"R^- Data collection run. 


2. The next digit 1s the grid number 


0- NO GRID 
F- GRID | 
2A GRID 2 
3- GRID 3 
4- GRID 4 


3. The next digit is the nose number. 
Nose 1 is at 0° roll angle. 


Each nose represents a 45° roll angle increment. 


4. The last two digits are the run number. 
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APPENDIX C. NUMERICAL DATA LISTING 


The following data was collected during the data collection runs. A complete list of 
all the runs and respective experimental conditions 1s found in the RUN MATRIX of 


Appendix B. 
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RUNFILE> WOA101 


AOA Cn Cy суп Сах Cpt Cri C 
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APPENDIX D. DATA ACQUISITION PROGRAM 


The following computer program was originally written by Sestak [Ref. 47] in 
BASIC and modified by the author for use with this experiment. This program follows 
a data aquisition shell program which allows the data acquisition program to control the 
Hewlett-Packard data acquisition hardware. The original program implemented auto- 
matic data input for angle of attack and tare values. Angle of attack is manually input 
to this program by the operator. 

Preliminary test runs indicated a larger than desired variation in the balance 
readings. The program was modified to average ten separate readings per channel and 
calculate the mean and standard deviation. Anv readings which were greater than one 
standard deviation from the mean were discarded. The program was also modified to 
allow the operator to recycle the balance reading routine without having to leave the 
program. This was necessary to enable the operator to reject all balance readings during 
a current read cvcle after viewing the data. 

Since the missile model pitched in the horizontal direction, only one tare reading 
was necessarv, eliminating the original routine for recording tare values at each angle 
of attack. Additional information for the program procedure is contained in the the 


program listing. 
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1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1150 
1140 
1150 
1160 
11υ 
1180 
1130 
1200 
1:10 
1220 
1250 
1240 
12510 
1200 
1270 
1280 
T2 
1300 
1919 
1-20 
1330 
1340 
1550 
το 
1570 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
15 
1540 
1550 


"Program to scan with the DMM and RELAY. MUX. 01 
‘This program was writen by T. SESTAK and modified by 
,P. RABANG for use with the VLSAM model. 


"This section after the SHELL program directs reading 

'the voltages from the balance, computes forces measured 

'by the strain guages, then stores the values in two arrays, 
"опе Ғог the TARE one for FORCE. This data file can then 
'be used for graphs or other displays. Each test run 

'Will generate a windtun.dat file which should be copied 
‘under another name before the next test run so that it 
will not be overwritten. 


‘dimension ‘arrays 
DIM READINGI[ 7] , FORCE[ 140,8] , TARE[ 8] , TREAD[ 7,10] 
COLOR 14,1,1 
CLS 


ΑΟΑΞΟ 

VALUE=5 

LOCATE 12,28 

PRINT'SETTING UP DATA FILES" 


"The program will write the data to several files. 


STATEFILES = "C: РСІВ WIND.HPC" 'stored in PCIB subdirectory 
DATAFILES 2 "C: WINDTUN. DAT" "stored on аі5с С 

DISKFILES = "A: WINDTUN. DAT" ‘stored on disc A 

BALANFILES = "C: MISSILE BALANCE. DAT "stored om disc С 
RELAY. SETTLING. TIME = .8 ‘800 ms 


LOCATE 16,35: PRINT'D ON E" 
CALL DELAY( VALUE) 


CLS: LOCATE 12,28: PRINT" INITIALIZING INSTRUMENTS" 
CALL INITIALIZE. SYSTEM( STATEFILES ) 

IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

CALL ENABLE. SYSTEM 

IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

τος ΝΡ 

CALL DELAY(VALUE) 


‘This part of the program is to preserve the data if 

'if the program is aborted in mid run. Parity errors 

'in the Hewlett Packard PC Instruments setup caused by 
‘electrical noise and undervoltage at NPS requires 

‘this. A voltage regulated, uninteruptible power source 
‘would ameliorate this problem. Just in case- this little 
'sequence allows reentry into the program and the data 
arrays with minimal inconvenience. 


CLS: LOCATE 12,20: INPUT'WERE YOU INTERRUPTED (Y OR N)"; A$ 
IF AS="Y" GOTO 1670 
1 


'The next three variables are counters in the arrays 
‘FORCE and TARE 


130 


1560 
1570 
1250 
1290 
1600 
1610 
1620 
1630 
1650 
1650 
1660 
1670 
1680 


1690 
1700 


1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1620 
1830 
1850 
1850 
1560 
18/0 
1880 
1890 
1900 
1910 
1920 
1950 
1940 


1950 
1960 
1970 
1980 


ТАТ 
ТКҮ = 0 


0 


| open the datafile so each scan can be recorded 
OPEN DATAFILES FOR OUTPUT ΔΒ 41 
CLOSE #1 
OPEN BALANFILES FOR OUTPUT AS #3 
CLOSE #3 
GOTO 1800 
OPEN DATAFILES FOR INPUT AS #1 
INPUT #1, TARE(1),TARE(2) ,TARE(3) , TARE( 4), TARE(5) , TARE(6), 
ТАКЕ 7) ;ТАКЕ(6) 
FOR X = 1 TO 140 
ΗΝ. FORCE(X,1) ,FORCECX ,2) ,FORCE(X,3) , FORCE(X,4), 
ВСК, 5). ВОВСЕСХ, 6), ЕОВСЕСХ, 7), РОВСЕСХ, 8) 
МЕХТ Х 
CLOSE #1 
GOTO 1800 
'AS is used as a marker for interrupted run sequences 
Potne progran, it is set to <>"'Y" so the 
uninterrupted Sequences are used unless otherwise directed 
A s= 
KEY OFF 


‘prompt to begin each scan or quit program if desired 
! 


685: LOCATE 12,10 
INPUT "TO START SCAN ENTER ANY KEY EXCEPT Q, Q TO QUIT"; ANSWERS 
IF ANSWERS = "Q" THEN GOTO 4730 


‘THIS ENTERS THE AOA FOR EACH TRIAL AND DISPLAYS IT ON THE CRT 


ms TOCATE 12,10 
PRINT "THE CURRENT ANGLE OF ATTACK IS "; AOA 


LOCATE 15,10: PRINT "INPUT THE ANGLE OF ATTACK (AOA) FOR THE 
NEXT TRIAL" 

INPUT AOA ' READING( 1) 

GOTO 1970 


КЕАрІМС( 1)=АОА 


1990 ' 


2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 


‘This variable is a marker in the iteration loop 
‘interaction equations for convergence. 
t 


CYCLE = 0 

‘This loop scans the pitch angle and 6 balance channels 

"and stores the values in the array READING 

‘Each channel is read ten times and averaged. 

‘The user may reject the current readings and input a new set. 


iol 


2090 
2100 
2110 


2120 
2120 


2140 


2150 
2160 
2170) 
2180 
2190 


! 


R tt 


Joey е 


CLS 
ο. пса ПТКЕСТ ВАГАМСЕ БКЕАрІМС5 «τσοὺ 
ος ος ος ον ος ος τις. 
PRINT" CHECK OF SYSTEM OPERATION 
PRINT 
" IN VOLTS N1 N2 S SZ A 
PRINT 
И дукикких 999599 N 99222929270 Jeezy 299727992 929776709276 
1 
"This file is for storing the direct voltage readings and averages. 
"The data file is continually appended. 
'The data is for further analysis of the direct voltage readings. 
OPEN BALANFILES FOR APPEND AS #3 


200059 


22 [0 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
23510 


2320 


29910 
2340 
23910 
2360 
2970 
2990 


2550 


2&00 


2&10 
2&20 
2&30 
2440 
2450 
2460 
2470 
2480 
2490 
2500 


FOR CNT = 1 TO 10 
FOR CHANNEL - 2 TO 7 


CALL OUTPUT(RELAY.MUX. 01, CHANNEL) 
IF PCIB. ERR <> 0 THEN ERROR PCIB. BASERR 


CALL DELAY(RELAY. SETTLING. TIME) 


IF PCIB. ERR <> 0 THENCERBRORCPONBPBSBASISEE 
CALL MEASURE( DMM. 01, READING[ CHANNEL] ) 
IF PCIB.ERR <> ο ΤΕΝ ΕΚΕΟΣ TO Ip pa P h 
TREAD(CHANNEL,CNT) - READING(CHANNEL) 


NEXT CHANNEL 
PRINT USING 
" +. HHEHABHE Ε.1ΗΗΗΗΗΙ - 


іНННННЬ +. ТТ 


+. НЕ 


+. EEN 


READING(2),READING(3),READING(&),READING(5),READING(6),READING(7) 


PRINT #3, USING 
"HHH + HHH 1. +. 


EE +. HHH 


+. НН 


+. HEHE; 


READING(1) ,READINGC2)  READING(3), READINCCA) | READINGCS) í D 


READING( 7) 
NEXT CNT 
! 


' CALL SUBROUTINE TO AVERAGE READINGS 


GOSUB 4790 
! 


PRINT USING 
"MEAN VALUE -.jJHHHHHE 1. +. 


НЕНННЬ +. HHHHH 


TOPIHHRHE 


tOBBHHHE'; 


READING( 2) , READING( 3) ,READING(4) ,READING(5) , READING(6) , READING(7) 


PRINT #3, USING 
"о НЫ + HHHHHE +. HHHHHI 


TOIHHHHHE +. HHHHHE 


+. НТ 


tHHHHHHE 


READING( 1) , READING( 2) ,READING( 3) ,READING( 4) ,READING( 5) , READING( 6) , 


READING(7) 
CLOSE #3 
PRINT" ": BEEP 


PRINT’ <CR> TO CONTINUE, "1" TO GET NEW READINGS" 


INPUT XYZ 
IF XYZ-1 GOTO 2040 


These equations take voltage readings from the balance, 
converts them to counts, then applys the primary force 
equations to the results. These values are applied to 


'the balance interaction equations. 
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Each channel has 


‘separate equations for positive and negative readings and 
2520 'тау have a "-" or "-" reading on any test run so the 
2530 "rather involved logic path below is my solution to the 
2540 ‘problem. For more information consult Calibration laboratory 
2560 guidelines ac NASA Ames Research Facility for TASK balances 
2350 
2530 UE dE AMET CONVERT SIGNAL TO FORCES "eme 


^ €: ! aub m as wap map ab m a m ap m apm arm ap m п 0 2a en ee 20 20 ee 20 20 ewe wes ese we se ee 20 seem, Бұрған тарта дала ла лалар ама ла жаа os ne alee we wn ee Кеө ө 
Ὁ 4Ν 4* «ν «ν 29 9% 6% 7% 7% 79 4% 6% 79 ОУ 4% 6% 4% 7% 79 79 ғ% 4% 2% 79 ғ% 4' 6% 79 7% ὅν 4 όν eee FE FC FD FED ED FE FD EX FB FH EY HY EL FR FR ED FD FH FY όν 
= 


2590 ' 

2600 'Direct balance readings are multiplied bv a scale factor 
2610 '5000000 then divided by the balance exitation voltage to 
2620 'get a readiing in COUNTS. The program will send each reading 
2630 'to the appropriate equation and convert to force or moment 
2660 then return to send the next reading for calculation 

2650 'The data acquisition system for using this program used an 
2660 ‘amplifier with 1000 gain. The scale factor is divided by 1000. 
2670 ' 

2680 УЕХ=5 ‘Excitation Voltage 
2690 Ni-READING(2)*5000! /VEX 

2700 N2-READING( 3)75000! /VEX 

2710 Si=READING(4)*5000! /VEX 

2720 S2=READING(5)*5000! /VEX 

2730 A=READING(6)*5000! /VEX 

2750 R=READING( 7)*416. 673/VEX 

2750 | 

2760 'send each reading to the appropriate equation 
2770 
2780 IF READING(2)»0 THEN GOTO 2870 ELSE GOTO 3020 
2790 IF READING(3)>0 THEN GOTO 2890 ELSE GOTO 3040 
2800 IF READING(“)>0 THEN GOTO 2930 ELSE GOTO 3080 
2810 IF READING(5)>0 THEN GOTO 2950 ELSE GOTO 3100 
2620 ТЕ READING(6)>0 THEN GOTO 2910 ELSE GOTO 3060 
2830 IF READING(7)>0 THEN GOTO 2970 ELSE GOTO 3120 


-- 


2850 ο ας POSITIVE FORMULAS 22922227272:76229292225797 


СО 1 = .0508615N1 - 5.4826E-09*(N1*NI) 
2960 GOTO 2/90 

29003082 = 7949211%Ң2 - 1./015Е-08%( 52642) 
2900 COTO 2800 

СО ША = . 01913097A - 7. 1962E-10*(A*4) 
8620 GOTO 2880 

О №. : - 03130951 - 3.81538-08*(Si*S1) 
2920 GOTO 2810 

ου = 17030366752 - 3. 86075-08" (52—52) 
2960 СОТО 2820 

2990 BR — .0030885*R + 2.5672E-09*(R*R) 
2980 GOTO 2190 


2990 
3000 κα: NEGATIVE FORMULAS м 
3010 ' 


ЗОВО ο = .05159J99NÀ ^ 1.7157E-08^(NT*N1) 
3030 GOTO 2790 

3040 EN2=. 047763*N2+8. 915299E-09*(N2*N2) 
3050 GOTO 2800 

30860 S= 0454529*A - 1. 3322E-09"( A*A) 


3070 
3080 
3090 
3100 
ο 
2120 
3130 
3140 
3150 
3160 
3170 
3180 


2190 
3200 


3210 
3220 


5250 
3240 
250 
3260 
3270 
3280 
3290 
3300 
3510 
5320 
3390 
3340 
350 
3360 
3370 


2560 
390 
3400 


3410 
3420 
3430 


3440 
3450 
3460 
3470 
3480 
3490 


3500 
5510 
η 


GOTO 2830 

ES1 - .032073*51 - 679316008 =-09-7 5st) 

GOTO 2810 ' 

ES2 = .031167%#S2 - 7.2517Е-09%(52%52) 

GOTO 2820 

ER = .0030908*R - 2.4769E-09*(R*R) 

1 

‚а heading for the iteration values 

PRINT" " 

PRINT seeded ο ος ος ος FORCE INTERACTION ITERATIONS 

Уза сат азал ыле 52 S 

PRINT" à CHECK FOR CONVERGENCE 

PRINT 

ο οποιο. №1 №2 51 S2 А RSS 
PRINT 

+ $ DEG POUNDS POUNDS POUNDS POUNDS  POUNDS  FT-LBS" 
PRINT 

H eee ЖУ: eee dee у зе γοῦν еее Sevededevede 
1 

'The loop that controls the balance interaction 

equations and allows a visual convergence check 

FOR I = 1 TO 10 


IF READING(2)^20 THEN СОТО 3270 ELSE 60103570 
IF READING(3)>0 THEN GOTO 3400 ELSE GOTO 3600 
IF READING(4)>0 THEN GOTO 3460 ELSE GOTO 3660 
IF READING(5)>0 THEN GOTO 3490 ELSE GOTO 3690 
IF READING(6)>0 THEN GOTO 3430 ELSE GOTO 3630 
IF READING( 7)>0 THEN СОТ Ә-Ә ЕББЕ ЕВ 2720 


КА κο ο. FORMULAS? 9 eR 

1 

XN12 EN1+. 0058036*N2+. 0041655*S1+. 058079*R-7. 1926E-07*(N2*N2 ) 
4. 0352Е-06%( 51%51)-.0006786%(КЭК) 

СОТО 3290 

1 


XN2= EN2+. 046218*N1-. 0028393*A-. 0081694*S1+. 0041463*S2+. 077279 
*R-6.8577E-07*(NI*N1)-1.7755E-05* CAXA) E2172 I9E-065( S 15:9 115505 220 
*(82*82)-.0019294*(R*R) 

GOTO 3300 


XA= EA+8. 6893E-04*N1+6. 0359E-04*S514+7. 7722E-05*52-. 11115*R 

46, 6537Е-07%(М1%М1)-4. 7936Е-06%(51%51)-671033Е-0650525520 0 pn 
*(R*R) 

GOTO 3330 


XS1 = ES1-6. 3459E-04*N1-. 11148*R+5. 535E-06*(N1*N1)+. 0024592*(R*R) 
GOTO 23310 


XS2= ES2-.0024237*N]14. 0022455* Act. 00667854 S To ο οτι ο κο 
*(N1*N1)4+1. 2072E-05*( A*A) -2. 7825E-06*(S1*S1)+, 0062217*(R*R) 
2009 5528 

1 


XR= ER-1.9928E-04*N2-2.5893E-04*8241. 1512E -07*(N2*N2) -5. 156E -08 
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3530 
3550 
3550 


%(52%62) 

GOTO 3730 

= οσο NEGATIVE FORMULAS ea S a sas oss sss 

XN1= EN1+. 010257*N2-. 0045396*S1-.04494*R+7. 9499E-07*(N2*N2 ) 


mle der h-O06*(S1%51)-.0003232"(R*R) 
GOTO 3290 


XN2= EN2+. 051778*N1-. 0044056*A-9. 038499E-03*S1-. 061125*R 
+5. 2897E-06*(N1*N1)41. 0467E-05%(A*A)-4. 8493E-07%(S1*S1)-. 0011773*(R*R) 
GOTO 3300 


XA- EA-.0021217*N1-*9. 1524E -04*N2-. 097148*R- 4. 2547E-06* (N1*N1) 
+4, 5846E-06*(N2*N2)-7. 5001E-04*(R*R) 
GOTO 2330 


НТ ЕСІ 007127551-. 0089235%А-. 05268%Қ-1. 2923Е-05%( Қ1%51) 
-4.0345E-05*(A*A)-9. 3969E-04*(R*R) 
5049075510 


ου οο το 0052619*N2-. 0072915*A*. 006856*S1-. 062581 
ОЮ 1TE-0720N15N1)-8. 6265E -06*(N2*N2) 3. 70545E-05*(A* A) T9. 983001E-06 
*(S1*$1)-*8. 0007E-OA* (R*R) 


GOTO 3320 

AR= ER+3. 5945E-04*N14+1. 5497E-07*(N1*N1) 
"Shift all the new variables back to the old name 
'for the next iteration 

N1=KN1 

N22XN2 

А=ХА 

S1=xS1 

S2=XS2 

ΕΞΣΕ 

! 

'A marker for the interations 

Cece = CYCLE + 1 


‘print the iterations to watch for convergence 
! 


PRINT USING 

i ## НЫЕ +I. +H d Ημ. tHE dE СЫННЫҢ ἩΗΗΟΊΗ; 
CYCLE ,AOA,N1,N2,$1,82,A,R 

NEXT 1 


INPUT "IF CONVERGENCE IS ADEQUATE ENTER Y, IF ANOTHER RUN 
IS DESIRED ENTER N'; ANSWRS 
IF ANSWRS = "N" THEN СОТО 3160 


NORMAL = N1 + N2 


SIDE = S1 + 82 

AXIAL = A 

PITCH = (N1-N2) * .0854 
ΚΟΠ, 5 (51-52) * .0698 
YAW = R*¥12 


525 


2990 
4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 


4210 
4220 
4230 


4240 
4250 
4260 


4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 


! 


TRIAL = TRIAL + Í 

INPUT; "IS THIS A TARE READING, Y OR N"; ANS 

IF ANS <> "Y' GOTO 4290 

COLOR το 1 

ТЕАТ S 

TRY = TRY + 1 

IF A$<>"Y" GOTO 4100 

CLS: LOCATE 12,10 

INPUT ENTER THE VALUE OF THE LAST TRY BEFORE INTERRUPT"; ITRY 


TRY = ITRY + 1 

TARE(1) - ТВҮ 

TARE( 2) = AOA 

TARE(3) = NORMAL 

TARE(4) - SIDE 

TARE(S) = AXIAL 

TARE(6) = PITCH 

TARE(7) = ROLL 

TARE(8) = YAW 

PRINT" " 

PRINT '*% % oe κκ Χα ΚΑΣ Δα TARE CALCULATION A k 

de ye de κ το oo 

PRINT 

ТРА AOA NORMAL SIDE ΛΑ PITCH ROLL YAW" 
PRINT 

Т j; DEG POUNDS POUNDS POUNDS ЕТ-185 FT-LBS  FT-LBS" 
PRINT 


11 4 


9%: И 


959292977: Toce 
'A loop to list all values so far 


FOR J = 1 TO TRY 

PRINT USING" HO HH HE — ÓHHIL1HE— ΠΗ ΗΝ ЗНН. ЗБЕ 
ІНЕ, Е Jik. HE"; TAREC1) , TAREC2) , TARE(3) , TARE(4) , TARE(5) , TARE( 6) , 
TARE( 7 ) , TARE( 8) 

NEXT J 


GOTO 4550 


“еее Я 95969672556 99299997 еее υσος 


ТЕ А$<>"У" GOTO 4340 
CLS ο ο σα. 
INPUT"ENTER THE ee OF THE LAST TRIAL BEFORE INTERRUPT"; ITRIAL 


TRIAL = ITRIAL + 
PORCH Gn = ΠΤΙ 
FORCE(TRIAL,2) = AOA 
FORCE(TRIAL,3) = NORMAL - TARE(3) 
FORCE(TRIAL,4) = SIDE - TARE(4) 
FORCE(TRIAL,S) = AXIAL - TARE(5) 
FORCE(TRIAL,6) = PITCH - TARE(6) 
FORCE(TRIAL,7) = ROLL - TARE(7) 
FORCE(TRIAL,8) = YAW - TARE(8) 


4420 ' 


4430 


‘print the values and store in file 


4440 ' 


4450 
4460 


4470 


PRINT" " 

PRINT''ve % ve ye ve ov x Ye X π FORCE CALCULATIONS 3 * * 
se ve κους X 

PRENDS 
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4480 
1 


+590 
tt 


&500 


ΙΑ, AOA NORMAL SIDE AXIAL PITCH ROLL ЖАУ» 


4510 ' 


4520 
4530 


4540 
4550 
4560 
4570 
4580 
4590 


4600 
4610 


4620 
4630 
4640 
4650 
4660 


4670 
4680 


4690 
4700 
4710 
4720 
4730 
8780 
8750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 


is DEG POUNDS . POUNDS POUNDS  FT-LBS  FT-LBS  FT-LBS" 
PRINT 
гие РКТ S ету 77272979796 Te n Tee a S ον) 
"а loop to list all values so far 
FOR J = 1 TO TRIAL 


PRINT USING" $$ — HH GS ННН НЬ СЫННЬЗНЬ tHE UE НН. ЗН) 
B oO duo. JS SFORGECJ, 1),FORCE(J,2) ,FORCE(J, 3), FORCE(J,4), 
Bu 509559 SEORCECJS69,;FORCE(J,7),FORCE( J,8) 

NEXT J 


‘Write the data to the output data files 
COLOR 14,1,1 
OPEN DATAFILES FOR OUTPUT AS f 
BSNESUSSFARECOIOSTARE(2);,TARE(3), TARE( 4) , TARE(C5) , TARECÓ), 
ERRE( 7), TARE(8) 
FOR X - 1 TO 140 

DONNEES I ЕСКСЕГХ, 1), FORCE(X,2) ,FORCE(X,3) ;FORCE(X,4), 
PORE S) SEORCE(X,6),FORCE(X,7),FORGCE(X,8) 
MEXT X 
CLOSE 51 


OPEN DISKFILES FOR OUTPUT AS 72 
ШЕШЕ 72- ТАБЕСІ), ТАКЕ(2),ТАКЕ(3),ТАКЕ(4) , TARE(5) , TARE(6) , 
ТАКЕ(7),ТАКЕ(8) 
FOR X - 1 TO 150 
ЕЕЕ О SPORGECX,1),FORCE(X,2),; FORCE(X, 3) ,FORCE(X,4), 
ВОВЕ, 5), РОБСЕ(Х, 6), РОКСЕСХ, 7) ,РОКСЕСХ,8) 
MEXT X 
CLOSE #2 


'Prompt for next scan 
INPUT "DO YOU WANT ANOTHER SCAN (Y OR N)'; ANSWS 
IF ANSWS <>"N" THEN GOTO 1790 


END 


'This subroutine averages the balance voltage readings 

'by computing the mean and standard deviation. 

‘Any readings less or greater than one standard deviation 
| 'are thrown out and a new mean is computed 


FOR CHANNEL Ξ 2 ΤΟ 7 


№=10 
SREAD=0 
SSDEV=0 
'Mean of balance voltage readings 
FOR CNT = 1 TO 10 


SREAD = SREAD + TREAD(CHANNEL,CNT) 
NEXT CNT 

MEAN = SREAD/N 

READING(CHANNEL) = MEAN 
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4950 
4960 
4970 
4980 
4990 
5000 
5010 
2020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5 bI 
321: 
ЭТ 0 
5140 


IFCN < 10) НЕМ бОТО БО) 


'Standard deviation routine 
ЕОК СУТ - 10180710 
DIF » TREAD(CHANNEL,CNT) - MEAN 
SDEV = DIF * DIF 
SSDEV = SSDEV + SDEV 
KEST CNT 
DEV - SQR(SSDEV/N) 
HI = MEAN TTF DEV 
LO = MEAN - DEV 
FOR CNT ου IO 10 
ARG = TREAD(CHANNEL,CNT) 
IF (ARG < HI) AND (ARG > LO) THEN GOTO 5110 
TREAD(CHANNEL,CNT) = 0! 
N =N =-= 1 
NEXT CNT 
NEXT CHANNEL 
RETURN 
END 
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APPENDIX E. COEFFICIENTS TRANSLATION PROGRAM 


This program was written by the author for use with this thesis experiment. The 


program reads a data file containing force and moment values and converts these into 
a coefficient form taking into account: 


Test conditions. 


Body configuration for blockage. 


Screen number. 


The results are written into a data file which are used to produce graphs and plots. 
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1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1155 
1200 
1210 
t220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1210 
1320 
1330 
1340 
1350 
1360 


1370 
1380 


1330 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 


1500 
1510 


' PROGRAM BY M.P. RABANG TO READ FORCE AND MOMENT VALUES FROM 
' & DATA FILE CREATED BY THE DATA ACQUISITION PROGRAM. 


COLOR 10,4,1 
DIM TARE[ 8], FORCE[140,8] , COEF[ 140, 9} 


CLS 
LOCATE 5 

INPUT" ENTER THE NAME OF THE INPUT FILE"; INFILES 
INPUT’ ENTER THE NAME OF THE OUTPUT FILE"; OUTFILES 
INPUT ENTER THE MINIMUM TEMPERATURE (F)''; TMIN 
INPUT"ENTER THE MAXIMUM TEMPERATURE (F)'; TMAX 
INPUT'ENTER THE BODY TYPE, A-*, В-0, С-Х "; BODS 
IF BODS-"A" THEN BODY-1: GOTO 1180 

IF BODS="'B'' THEN BODY=3: GOTO 1180 

IF BODS='"'C'' THEN BODY=2: GOTO 1180 

GOTO 1130 

INPUT"ENTER THE SCREEN NUMBER, 0-4';SCR 

IF SCR=0 THEN GOTO 1250 

ΙΕ SCR=1 THEN 0-215995 C ROI? CU 

IF SCR=2 THEN Q=14. 882: GOTO 1260 

IF SCR=3 THEN Q=16.38: GOTO 1260 

IF SCR=4 THEN Q=15.605: GOTO 1260 

GOTO 1180 

INPUT"ENTER THE WIND TUNNEL VELOCITY (cm H20)"; PH20 
INPUT"DO YOU WANT TO VIEW OUTPUT"; PANSS 
INPUT"DO YOU WANT A HARDCOPY"; ANSS 


' OPEN THE DATA FILE SO EACH SCAN IS RECORDED 
TRANSFILES="C: MISSILE TRANS. DAT" 

OPEN OUTFILES FOR OUTPUT AS #2 

CLOSE #2 


' READ THE FORCE VALUES FROM THE DATA FILE 

OPEN INFILES FOR INPUT AS #1 

INPUT #1, TARE(1),TARE(2) ,TARE(3),TARE(4),TARE(5) ,TARE(6), 
TARE( 7) , TARE(8) 

FOR X = 1 TO 140 

INPUT #1, FORCE(X,1),FORCE(X,2),FORCE(X,3),FORCE(X,4), 
FORCE(X,5),FORCE(X,6) ,FORCE(X,7) ,FORCE(X,8) 


IF FORCE(X,1)20! THEN GOTO 1410 
NEAT X 

CLOSE ІН 

m PANSS<>"Y'' THEN GOTO 1610 

' DISPLAY FORCE VALUES 

COLOR 14,1,1 

CLS 

PRINT FILE NAME: '; INFILES 

PRINT" " 

PRINT x voveo ve veo we ve de de de ve ve ve X FORCE READINGS 
еме с se Se se Se ος Se u 22 0 20 50 
PRINT" " 

PRINT 


"TRI Ali eos NORMAL SIDE AXIAL PITCH 


140 


1520 


1530 


1540 
1550 
1560 
5 70 


1580 


ROLL YAW" 


PRINT 

um ył DEG POUNDS POUNDS POUNDS FT-LBS 
FT-LBS FT-LBS" 

PRINT 

ο. r i cas ом 5 2v ἐν όν 28 TX 2v ΘΝ ον ον ἐν ἐν αν NL ON CRON ON ON ON ON πο ος οὐ ος ον οὐ οἱ ον 
1 

FOR J = 1 TO 140 


IF FORCE(J,1)=0 THEN GOTO 1590 

PRINT USING" дай + HEH HE HHE HE HEEL HH 
О о НЫ: — HH HHRHE'; FORCECJ,1) ,FORCECJ,2) ,FORCE(J,3), 
F@Rer(),4),FORCE(J,5),FORCE(J,6), FORCE(J, 7) ,FORCE(J,8) 

NEXT J 


1590 ' 


1600 
1610 
1620 
1630 
1650 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1790 


INPUT'ENTER <CR> TO CONTINUE"; INPTS 


' BEGIN COEFFICIENTS CALCULATION 
A=. 0167 

Me=3.719E-07 

TAV=( TMIN+THAX) /2! 

RHO=2116. 22/(1545*(459, 7+TAV) ) 
IF SCR>O GOTO 1740 

ШЕ 2. 040-720) /(. 93*RHO) 
VEL=SQR( VE) 

RED=( RHO*VEL* (1. 75/12) )/MU 

Q=( RHO*VEL* VEL) /2! 

GOTO 1800 


VEZ2*Q/RHO 

VEL=SQR( VE) 
RED=(RHO*VEL*(1.75/12))/MU 
PH20=Q*. 93/2. 046 


~ WRITE ΜΕ COEFFICIENTS TO THE OUTPUT FILE 
OPeM OUIFILES FOR OUTPUT AS jf 

ОРЕ TRANSFILES FOR APPEND AS #3 

Ur UNE 5:35 - OUTFILES 

FOR X = 1 TO 150 

poser - 2 TO 8 

IF FORCE(X,1)=0 THEN GOTO 2040 

' ROUTINE TO CORRECT THE DYNAMIC PRESSURE FOR BLOCKAGE 
' AND BODY TYPE 

COEF@. ,1)=FORCE(X,1) 

COEF(X,2)-2FORCE(X, 2) 

ALPHA - ҒОКСЕ(Х,2) 

IF FORCE(X,2) < 0 THEN ALPHA - ABS(FORCE(X,2)) 
IF FORCE(X,2) > 90 THEN ALPHA = 180-FORCE(X,2) 
IF BODY=1 THEN EPS=. 0000126*ALPHA+. 007759 

IF BODY=2 THEN EPS=. 0000101**ALPHA+. 007759 

IF BODY=3 THEN EPS=. 0000908*ALPHA+. 007759 

D A*Q*( 1+( 2*EPS)) 

COEF(X,Y) = FORCE(X,Y)/D 
COEPCX,9)=COEF(X,4)/COEF(X, 3) 

NEXT Y 
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IF FORCE(X,1)=0 THEN GOTO 2040 

WRITE #2, COEF(X,1),COBE(X 2) COEE(CX 3) COEF( X О Wa pm КОО n 
COEF(X,5) ,COEFCK, 6) Cen GOD S 

WRITE #3, COEF(X,1),COEF(X,2) ,COEF(X,3) ,COEF(X,4) ,COEF(X,9), 
COEF(X,5),COLFCX,6) ,CORRGX ο ος ο 

NEXT X 

CLOSE #2 

CLOSE 3 

IF PANSS«»"Y" THEN GOTO 2300 


' DISPLAY ROUTINE 
COLOR MI 


CLS 

PRINT'FILE NAME: "; OUTFILES 

PRINT" " 

PRINT'AVERAGE TEMPERATURE (F) "; ТАУ 
PRINT'WIND TUNNNEL VELOCITY (FT/SEC) Pe VEL 

PRINT WIND TUNNNEL VELOCITY (cmH20) "; РН20 
PRINT'AIR DENSITY (LBM/FT3) И; ВНО 
PRINT'REYNOLDS NUMBER "; ВЕБ 
PRINT'DYNAMIC PRESSURE (LB/FT2) ο 

PRINT" " 

PRINT'* se ok c % e (ж FORCE COEFFICIENTS 
η... ο”. e 

PRINT" " 

PRINT 

"TRIAL ο NORMAL SIDE ASTAL PITCH 
ROLL YAW" 

PRINT 

ΠΤ ος ους тее sevesecedevevieve еа a i a aa еее ое ΛΕΩ 
το οσο ολ 

ΕΡΕ X = 1 TO 140 


IF COEF(X,3)=0 THEN GOTO 2280 

PRINT USING" ΗΕ HH. db +##. ff RH THHHE НЕ. ННІ 

THO. BERE .Η ΗΠ οὐεε; ΠΕ, ЕНЕС ЕС я 
COEF(X,4),COEF(X,5),COEF(X,6) ,COEF(X,7), COEF(X, 8) 

NEXT X 


PRINT" S 

IF ANS$<>"Y" THEN GOTO 2640 

' HARDCOPY ROUTINE 

LPRINT' INPUT FILENAME: ''; INFILES 


LPRINT OUTPUT FILENAME: "; OUTFILES 

LPRINT" " 

LPRINT INITIAL TEMPERATURE (F) ". TMIN 
LPRINT'FINAL TEMPERATURE (F) "; ТМАХ 
LPRINT'AVERAGE TEMPERATURE (F) ΝΡ 
LPRINT'WIND TUNNEL VELOCITY (FT/SEC) Ευ 
LPRINT'WIND TUNNEL VELOCITY (cmH20) "; РН20 
LPRINT AIR DENSITY (LBM/FT3) "5 RHO 
LPRINT'' REYNOLDS NUMBER "; RED 
LPRINT'' DYNAMIC PRESSURE (LB/FT2) "0 
LPRINT S 

LPRINT" " 

LPRINT ese e eue ο ο... б FORCE READINGS 


2490 


2500 
2510 
2520 


2530 
255 

2550 
2560 


2570 
2580 


23590 


2600 
2019 
2620 


2630 
2640 
2650 
2660 
2670 


ә yo se ce ve se te ss € se se ve ves 95 ox 


LURINT" Ἡ 


LPRINT | 
" TRIAL AOA NORMAL SIDE AXIAL PITCH 
ROLL YAW" 
LPRINT 
в DEG POUNDS POUNDS POUNDS FT-LBS 
ου. PIERS 
Teele 

тай I ο ἄν ος κος ος ο ο σος ος eye eee eee 959292679277 
WOOD as την 


FOR J » 1 TO 140 

IF FORCE(J, 13-0 THEN GOTO 2540 

μμ USING Geese ΕΓ. η. +G. HHH  +##. Hi 

НР. И THE. НЫН ο = FORGE( J.,1),FORCE(), 2) ,FORCE(J, 3), 
FORCE), ESO CECI 5), FORCE(CJS6),EORCECJ 7) EORCE(J ,8) 

NEXT J 

БРАТ" 

Leer 

LPRINT'* д: 5-5: Be Ws ae as We 9 X FORCE COEFFICIENTS 


LPRINT" ' E 

LPRINT 

Pore VAL AOA NORMAL SIDE AXIAL PITCH 
ROLL YAW" 

LPRINT 

По узок Ут и Е Ἐεοετεητος ΤΣ еее 


о κιν. G * г!! 
22:727::%272%777 ыса eae жән v 


ШЫК Ж = 1 ТО 150 

IF COEF(X,3)=0! THEN GOTO 2640 

LPRINT USING" Ны. Н THES4HHH) OO НН. ННН 
TÉ. 903 HH. i НН τι. 34H12": COEFCX,1), COEF(X,2) , COEF(X,3), 
CRER. 4), COEF(X,5) ,COEF(X,6), COEF(X, 7),COEF(X,8) 

NEXT X 


INPUT'IDO YOU WANT TO EXIT THE PROGRAM"; AANSS 


IF AANSS<>"Y" THEN GOTO 1060 
END 
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